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LIST OF SYMBOLS AND ABBREVIATIONS 



Amn 


solar insolation absorbed in the layer (m,n) 


a (m,n) 


Manabe-Moller absorptivity function 


ABA 


absorptivity of the troposphere 


ABG 


fractional absorptivity of solar insolation by 
earth* s surface 


ALB 


earth-atmosphere system albedo 


ATRAN 


transmissivity of the troposphere 




Stefan-Boltzmann blackbody flux at T^ 


BALB 


24-hour averaged radiational balance at earth's 
surface 


BALk 1 k 2 


24-hour averaged radiational balance for layer 
(k l' k 2 } 


BALT 


24-hour averaged radiational balance at tropopause 


C 


carbon dioxide layer absorber mass 



cal cm min calories per centimeter squared per minute 



CL 


total opaque cloud cover 


CL i 


fractional cloud amount for layer: 
I = 1 in 600 to 400 mb; 

I = 2 in 900 to 800 mb 


E 


East longitude 


e 

X 


vapor pressure at top of constant flux layer 


F (A) 


solar insolation subject to water vapor absorption 
only 


F (2) 


effective solar insolation at tropopause 


FADJ 


total incoming insolation at top of atmosphere 


Fk l k 2 


net infrared flux divergence in layer (k^ f k 2 ) 
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V 


net infrared flux at level k 


FNWC 


Fleet Numerical Weather Central 


F (S) 


solar insolation subject to Rayleigh scattering only 


f t 


net IR flux to space 


g 


2 

gravity = 9.8067 m sec 


f 


multiplicative factor for tuning cloud reflectances 


h 


hour angle 


H 


height of homogeneous atmosphere; 24-hour averaged 
hour angle 


HL 


population of gridpoints north of 25N 


I 


abscissa grid location 


IA10(m,n) 


solar insolation absorbed at surface with cloud 
condition (m,n) 


IS10(m,n) 


solar insolation at surface subject to Rayleigh 
scatter with cloud condition (m,n) 


J 


ordinate grid location 


k 


pressure level used in this study equal to 10a 


i • - 1 

ly min 


langleys per minute 


N 


North latitude 


P 


pressure 


p k 


pressure in millibars (mb) at level k 


q k 


mixing ratio at level k 


Qk l k 2 


24-hour averaged solar warming in layer 


QAVE 


24-hour averaged insolation at the tropopause 


2 N 


solar net insolation at level k = 0 


r 


Bowen ratio; actual earth-sun distance 


r 

m 


mean earth-sun distance 
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R 


net radiation balance at the surface 


R 

a 


mean radiative cooling rate in troposphere 


R d 


universal gas constant 


REF 


total insolation reflected back to space 


REFA 


F(A) insolation reflected back to space 


REFS 


F (S) insolation reflected back to space 


RH 


relative humidity 


R t 


mean radiative energy gain (loss) rate at ocean- 
troposphere system 


r n 


radiative net flux at level k = 0 


s 


South latitude; effective solar constant 


s 

o 


heat storage in oceanic water mass 


S 

a 


heat storage term for the troposphere 


Study A 


Spaeth's Thesis; using winter data (see references) 


Study B 


Meyers' Thesis; using spring data (see references) 


Study C 


Beahan's Thesis; using summer data (see references) 


Study D 


Warner's Thesis; using autumn data (see references) 


T k 


temperature at level k 


TR 


population of gridpoints 20S-to-25N inclusive 


TRAN 


total insolation incident at the earth's surface 


T 

X 


temperature at the top of constant flux layer 


U 


water-vapor layer absorber mass 


W 


West latitude 


W (m,n) 


cloud fractional weight for cloud condition (m,n) 


Z 


Zenith angle 


a (G) 


surface albedo 
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a(R) 

6 

A (ALB) 
e 

WC 



A 

7T 

P 

a 






Rayleigh clear sky albedo 

solar declination angle 

difference between ALBMOD and ALBRAS 

emissivity due to water and carbon dioxide absorber 
mass at indicated layer 

potential temperature at level k 

longitude 

surface pressure? pi = 3.1416 
density 

sigma pressure level used by FNWC, normalized 
to surface pressure 

latitude 
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I. INTRODUCTION 



This thesis is a refinement of previous radiational models described 
by (A) Spaeth (1975), (B) Meyers (1975), (C) Beahan (1975) and (D) 

Warner (1974) for use in the Fleet Numerical Weather Central (FNWC) pre- 
diction system. This study has as a primary objective, the comprehensive 
re-examination of the radiational physics in layers comprising the ocean- 
atmosphere system. FNWC atmospheric soundings defined at constant pres- 
sure levels for the four mid-seasonal dates of the "data year" 1973-74 
as previously examined in studies (A,B,C,D) were utilized in this study. 

The application of the FNWC radiative model may be made to any scale 
of analysis for which there is adequate resolution of the temperature 
and moisture data in the vertical. In the horizontal, the reliability 
of the data used here is consistent with that of the FNWC interpolation 
to gridpoints in the analysis procedure. Temperature and dew-point data 
in radiative soundings are typically reported to the nearest tenth of a 
degree. The radiative computations made here are applied to FNWC grid- 
points and are designed to make a one-hour forward-time step at gridpoints 
in the FNWC primitive equation forecast model, with special adaptions to 
0-levels . 

The specification of cloud amounts in two designated layers, one at 
a mid-level and the other at a low-level, has an important influence on 
the radiative-model dispositions, both in the short- and long -wave 
spectral regions. Initially in (A,B,C,D), the specification of the 
fractional amounts of CL^ and CL^ had been based on large-scale formula- 
tions developed by Smagorinsky (1960) , but during the course of these 
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studies it was found more realistic to modify initial CL-values to 
CL* = 2/3 CL. This reduction in cloud amount was also used here as it 
had been chosen to prevent CL* from exceeding unity and to afford 
better agreement with satellite climatology, such as planetary albedo , 
for compatible data periods . 

With the reduced cloud coverages CL^ ' and CL^ 1 from the earlier 
studies (A,B,C,D), it was possible to obtain reasonably close agreement 
in the computed terrestrial net flux at the top of the atmosphere and 
that observed for comparable NIMBUS III subsatellite points and data 
periods (Raschke, Von der Haar, Bandeen and Pasternak, 1973). However 
even with the reduced cloud coverages CL^ 1 and the computed 

planetary albedo remained generally excessive, particularly in tropical 
latitudes. Hence it became a major objective of this particular study 
to modify empirically the reflective capability of the cloud layers. 

This was done by defining a general factor f so that the initial choices 
(after Rodgers, 1967) of cloud reflectances R were modified to R* where 

R* = fR o 

Systematic substitution of the cloud reflectances R, wherever they 
entered the solar disposition equations, by R* then led to a relation- 
ship between the global albedo and f. Utilization of the least squares 
technique to minimize the differences between satellite and model albedos 
over a geographic sample of points led to best-fit value of the "tuning- 
factor" f. Separate values of f were deduced by least squares for each 
season, and subselections were deduced for the tropical and extra tropical 
areas , respectively . 
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The modified solar cloud-reflectances improved the agreement 
between satellite and radiative model albedos in both geographic 
areas insofar as net incoming insolation was concerned. The terres- 
trial net flux at the top of the atmosphere was unaffected by the 
choice of f , while the use of CL^ ' and CL ^ 1 as specified gave good 
agreement with satellite terrestrial net flux data over the geo- 
graphic range and the time-scales concerned. 



15 



II. DATA PREPARATION 



A. DATA FIELDS 

1. General Considerations 

The initial temperature and humidity data used in this study 
were arranged in the form of soundings taken along four oceanic meri- 
dians (Fig. 1) of the Fleet Numerical Weather Central (FNWC) Northern 
Hemisphere mid-seasonal analyses for 16 October 1973, 16 January 1974, 
16 April 1974 and 16 July 1974. Oceanic locations for these computa- 
tions were chosen because: 

p 

(a) constant a-surfaces (where 0 = — ) of the FNWC primitive 
equation system are nearly identical constant pressure levels. 

(b) The maritime-area soundings are more likely to be system- 
atically representative of the set of zonally-distributed gridpoints 
than over land. 

The three meridians selected over the Pacific Ocean were lo- 
cated at 125W (25 soundings), 170W (25 soundings) and 145E (17 sound- 
ings) . The Atlantic Ocean meridian was located at 35W (26 soundings) . 
This method of selecting "soundings" along the indicated meridians 
of the FNWC polar stereographic map made it unnecessary to employ spa- 
tial interpolation between original data gridpoints along the meri- 
dians. Data along line 3 in the Pacific was not extended southward of 
gridpoint (9,55) because they fell over land masses (New Guinea and 
Northern Australia) where the surface temperatures and other sounding 
features were unrepresentative of the oceanic values. 
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'Figure 1. FNWC polar stereographic grid and meridians (lines 1, 2 , 

3 / and 4) selected for study. The longitudes A are shown 
for each meridian as well as the extent considered of each 
' meridian. 
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2. Data Treatment 



a. Original Soundings and Modifications 

The gridpoint soundings were taken from the original FNWC 
63-by-63 Northern Hemisphere analyses of T(p) and of T - T^ (the dew 
point depression) at standard pressure levels up to and including 
p ^ 100 mb. Examples of such original soundings were shown as Table 
1(a) in Meyers (1975, p. 24). Subsequently each original FNWC sound- 
ing was transformed — in previous studies (A, B, C, D) of this series 
— into what has been termed the radiative sounding having the format 
shown in Table I. The data levels of the radiative sounding contains 
essentially the five FNWC predictive O-levels (dotted levels in Fig. 2) . 

At each gridpoint selected, the original FNWC humidity 
soundings were given in the form of five dew point depressions over 
the analyses levels from 925 mb to 400 mb. At the surface (level 
k = 10) the standard instrument level vapor pressure, e a ^ r > was trans- 
formed into the surface mixing ratio, q^, by means of 

q_ n = 621.97 (e . /1000) . (2-1) 

lu ai r 

To obtain radiative soundings as in Table I, it is neces- 
sary to have water vapor and CC> 2 absorber masses at certain required 
k-level boundaries (Fig 0 2) . All radiative soundings in this study 
start at sea level with the approximation of surface pressure 
7T = 1000 mb. Therefore, the eleven k-levels correspond closely to the 
FNWC levels P = 1000., 900., 800., ...., 200., 100., 0.0 mb and in 
turn to C k = 1.0, 0.9, ...., 0.1, 0.0. 
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TABLE I. Example of a radiative sounding at gridpoint (1,1) for 
16 April 1974 with mixing ratio listed at odd k-levels (Fig. 2) . 
Additionally, water-vapor and CO- absorber masses are also listed 
as these parameters have been modeled in the radiative theory pre- 
sented in this study. 



Pressure 

(mb) 


Temp 

(°C) 


Mixing 

Ratio 

(gAg) 


Absorber 
Water Vapor 

(gm/cra^) 


Masses 
CO 2 

(cm/ cm 


1000 


25.60 


17.10 






900 


16.09 


12.08 






800 


11.21 




2.26 


45.53 


700 


5.00 


6.21 






600 


-2.10 




3.23 


83.53 


500 


-11.50 


2.57 






400 


-23.20 


1.35 


3.55 


113.35 


300 


-38.00 


0.65 






200 


-56.60 




3.61 


133.99 


100 


-80.70 


0.03 


3.61 


138.67 


0 


-80.70 




3.61 


143.35 
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I 



Lev-sl- 

Idcniifers 



U(0), C(0) 



O- level 

0=0 



Layer (0,2) 1* 



T,q- 



- 0.1 



0(2) , C(2 ) 



0 = 0.2 






10 



U(10), C(10) 



a=1.0 



Figure 2. Five-layer radiative sounding used in this study. Levels 
are identified by their values on the k-scale, while layers are identi- 
fied by their level boundary indices in parentheses, e.g. (8,10), 
Pressure-scaled water vapor and CO^ mass increments Au and AC, respec- 
tively are integrated relative to the surface and the resulting U and C 
are carried at even levels. The temperature T is retained at all levels. 
Amounts of clouds CLj and CL2 in the layers shown have been para- 
meterized for consideration of their radiative effects. 
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b. Radiative Temperature Profiles 

The gridpoint temperatures were listed at each mandatory 
level of Table I between 1000.,..., 100 mb (i.e., between k = 10,..., 
1) . The temperature was assumed to be isothermal from 100 mb to 
0.0 mb. The temperature at level k = 10 was set equal to the FNWC 
listed sea-surface temperature. The radiative sounding temperatures 
for the remaining k-levels were obtained from either their correspond- 
ing listed temperature-level or by a three-point Lagrangian interpola- 
tion scheme [Eq. (2-1), Spaeth, 1975], to level k when the listed 
FNWC temperature profile did not include the value T^. 

c. Radiative Moisture Profiles 

Similarly, the moisture profiles of Table I have been ob- 
tained by an interpolative procedure over the original-level FNWC mix- 
ing ratios to those required at k-levels in a manner analogous to that 
discussed by Spaeth (1975, pp. 29-31) . 

d. Pressure-Scaled Absorber Masses 

The pressure-scaled water vapor and the carbon-dioxide 
scaled absorber masses were calculated for the six even numbered k- 
levels (Fig. 2) using Eqs. (2-8, 2-9, 2-10 and 2-11) respectively as 
outlined by Spaeth (1975, pp. 31-33) . These equations use essentially 
the mixing ratios of water vapor and of CC> 2 at odd k-levels. 

B. CLOUD PARAMETERIZATION 

The relative humidities (RH) at levels k = 5 and k = 9 are used 
in the calculations of the fractional cloud covers CL^ and CL^ in 
layers (4,6) and (8,9) respectively (Fig. 2). The equations for 
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parameterization of the two fractional cloud amounts are as follows: 



CL = 2/3 [2.0 (RH ( 5) ) - 0.7] 



( 2-2a) 



CL 2 = 2/3 [3.33 (RH (9) ) - 2.0] 



(2-2b) 



The bracketed part of the equations (after Smagorinsky, 1960) were 
reduced by the 1 /3 factor in an attempt to tune cloud amounts to obtain 
albedo values in closer agreement with the recent satellite radiational 
climatology of Raschke et al., (1973). The "2/3 CL" parameterization 
set forth in Eq. (2-2) is used here for estimating large scale radia- 
tional effects only. Thus small-scale convective activity and ci priori 
climatological effects were not considered in specifying the form of 
Eq. (2-2) . Tuning of the model-albedo values of this study by varying 
cloud reflectance coefficients with respect to season and latitude will 
be discussed in Section IV. 

C. CLOUD-AREA COVERAGES 

Fractional cloud amounts, CL^ and CL^, were computed at each grid- 
point for levels k = 5 and k = 9 by Eqs . (2-2a,b) . In addition the 
gridpoint area may be thought of as broken into random fractional seg- 
ments of size 



wherein there is a combination of clear-over-clear segments in the 
layers. Similarly, the gridpoint area has the fractional area of cloud 
coverage 



W (0,0) = (1-CL 1 ) (1-CL 2 ) 



(2-3a) 



W(l,l) = CL 1 *CL 2 



(2- 3b) 
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of an upper overcast amount overlying a lower overcast. Likewise the 
area-combinations of overcast over clear, and clear over overcast 
areas in two layers, may be visualized as occurring with the weights 

W (1 , 0) = CL *(1-CL 2 ) (2-3c) 

and 

W(0,1) = (1-CL ) *CL 2 (2-3d) 

respectively . 

For radiational computations it was useful to carry the relative 
weights or fractions of the gridpoint area exposed to the specified 
cloud- layer combinations. Henceforth, the symbols denoted by W(0,0), 

W (1 , 1) , W(1,0) or W(0,1) indicate the fractionally overcast (1) or 
clear (0) cloud-area combinations in the indicated layers (Fig. 2) , 
with the first index 1 or 0 referring to layer CL n , k = 5, and the 
second to CL^, k = 9. 

The usefulness of the cloud-area weighting device will be clarified 
in Sections III and IV, where the procedures for the terrestrial and 
solar radiational computations are discussed and the results are 
summarized over the set of soundings. 

A measure of the effective cloud-cover area which has been found 
useful in previous radiational studies has been the total opaque cloud 
cover , CL, referring to the amount of thick cloud cover overhead re- 
gardless of the level. For the cloud model presented here CL may be 
expressed as 

CL = CL 1 + CL 2 - CL *CL 2 . (2-4) 
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III. TERRESTRIAL RADIATION 



A. THEORETICAL AND EMPIRICAL BASIS 

Empirical formulas were developed by Sasamori (1968) for flux 
emissivities in the atmosphere associated with computations for the 
radiative balance requirements of the NCAR General Circulation Model. 
Sasamori derived the empirical emissivity formulas for v/ater vapor 
and CO^ by comparison with the theoretical values built into the 
Yamamoto Radiation Chart (1952) . The Yamamoto chart has proved to 
be quite accurate for numerical checks of the Sasamori emissivities. 
This chart was also used in the previous studies (A, B , C, D) as a 
systematic guide for integration of the radiative transfer formulas 
developed by Martin (1972, 1975) , who adapted the Sasamori emissiv- 
ity formulas to the particular layers of interest in the gridpoint 
computations of the FNWC primitive equation model (Fig. 2) . 

The essential long-wave (IR) net-flux parameters required for 

use in this study are the following: 

* 

F = IR net flux at earth, k = 10 

* 

F = IR net flux at level k = 8 

o 

* 

F = IR net flux at level k = 6 
6 

★ 

F. = IR net flux at level k = 4 
4 

* 

F^ = IR net flux at level k = 2 
In addition the IR net-flux divergence coolings to be computed at 
each gridpoint are 

F810 = IR net-flux divergence in the layer (8,10) 

F68 = IR net-flux divergence in the layer (6,8) 
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F46 = IR net-flux divergence in the layer (4,6) 

F24 = IR net-flux divergence in the layer (2,4) . 

In the Radiation Balance Studies in the series A, B, C, D only 
F610 and F26 were computed because time restraints in the present 
FNWC operational heating package have prevented the use of greater 
resolution in the vertical. Here, four flux divergences are computed 
in order to examine more closely the variability of the flux diver- 
gences over the layer thicknesses reduced to approximately 200 mb 

each. To compute these four flux divergences it was necessary to 

* * 

utilize additional formulas for F and F as developed by Martin 

o 4 

(1975) . 

In order to make IR net-flux calculations along the path of in- 
tegration, there must be a physically sound representation of the 

emissivity (£ ) as a function of both water vapor and C0_ absorber 

wc 2 

masses in layers along the sounding. For a complete discussion of 
the emissivity formulas used in the quadrature scheme, refer to 
Spaeth* s Appendix A (1975) . 

B. NET FLUX FORMULATIONS 

The radiative sounding as depicted in Table I was computed as 

the combination of parameters U(k,10), C(k,10) and T for each re- 

K 

quired level, k = 10, 8, ... 1,0. Cloud parameters CL^ and CL 2 were 

also computed by Eq. (2-2) at each gridpoint and in general are both 

non-zero. The grid area was then considered to be composed of areal 

fractions (weights) defined in Eqs. (2-3a,b,c,d) and denoted by the 

symbols W(0,0) , W(l,l), W(1,0), W(0,1). 

* 

The composite net flux F^ Q (CL^, CL^) at level k = 10 at each 
gridpoint is then constructed by using the appropriate weight factors 
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•k 

to multiply the reference net flux F^ computations defined for the 
four special cloud-cover cases already defined in Section II. C: 

F 10* (0 ' 0) ' F 10* (1 ' 0) ' F 10* (0 ' 1) ' F 10* (1 ' 1) 

It therefore follows that 



F 10* (CL 1' CL 2 ) = W (°' 0 ) F 10 * (0 ' 0) + W ( 1 '°) F 10 * (1 ' 0) 

k k 

+ W(0,1)F 10 (0,1) + W(1,1)F (1,1). (3-1) 

★ 

The reference net fluxes of Eq. (3-1) are associated with (1) 

clear skies in both layers, (2) overcast in the upper layer only, 

(3) overcast in the lower layer only and (4) overcast in both layers 
respectively. 

Spaeth (1975) has listed these reference net flux formulations 
in his Eqs . (3-6), (3-7) and (3-8). Using the definitions of W(0,0), 

W(1,0), W(0,1) and W(l,l), F^* (CL^, CL^) Can s ^ own to assume the 
form 



F 10* <CL 1' CL 2> = [1 - CL 2 ^ ^ F wc ^ 8 ' 1C1 ^ ^ B 10 _B 8^ 

+ (e wc (8,10) + E wc (6,10))(B 8 -B 6 )]} 

+ (l-CLo ) (1 -CLh ){b -.5[(£ (6,10) 

+ £ wc ( 4 ' 10))( V B 4 ) + (£ wc (4 ' 10) 

+ £ wc ( 2,10))(B 4 -B 2 ) + (E wc (2,10) 

+ e wc (l,0))(B 2 - Bl ) + E wc ((0,10) / T 1 )*B 1 ]} 

+ CL 9 { (B — B ) [1.5c (9,10)]} . (3-2) 

10 9 WC 

Here 



B 



k 



1.170403 x 10 7 T 4 
k 



(3-3) 
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is the Stefan-Boltzmann blackbody flux in langlies per day. 

Further, £ (U, ,C . ,10) is the combined water-vapor and CCt emissiv- 

wc k k 2 

ity along the path from level 10 to level k. This emissivity is con- 
sidered by Sasamori to be temperature independent for T >_ 210K, whereas 

£ represents the temperature dependent emissivity applicable for 
wc 

T < 210K. [See pp. 136-137, Spaeth (1975); Sasamori (1968)]. 

* 



The formulas for F (CL , CL ) , k = 2,4,6, and 8 have been developed 

.K _L 2 

by Martin (1975) in a manner analogous to the derivation of the weighted 
* 

F . The results are reproduced as the following equations: 



F„ = 



[1 


- CL 1 


] 


{B 8 


-.5[e (6 

wc 


,8) 


( V 


-v 


+ 


<£ wc 


(6 


,8) 


+ £ 

wc 


(4, 


8)) 


( V 


-v 


+ 


(£ 


(4 


,8) 


+ £ 


(2, 


8)) 


(B.- 


-B„) 




wc 


WC 


4 


2 


+ 


(£ 

wc 


(2 


,8) 


+ £ 

WC 


(1, 


8)) 


(B 2- 


"V + 



v/c 



r l 



+ CL ± [l-.5e wc (6,8)] (B Q -B 6 )+ CL X (l-CL^ [ 1- . 5^ (8, 10)] (B^-Bq) 



+ (1-CL.) (1-CL 0 ) [1-.5S (8,10)] (B n -B_) 

1 2 wc 10 8 



(3-4) 



F 6 * - tl-CLj) (B 8 -.5[e uc (6,8)(B 8 -B 6 ) + (4, 6) (B^) 

+ <E WC <4 ' 6 > + £ wc <2 ' 6)) ( V B 2> + < E „c <2 ' 6 > 

+ £ wc (l,6)) (B 2 -B 1 ) + £ wc ((0,6) ,T 1 )*B 1 ]} 

+ (1-CL. ) (1-CL.) {(B in -B Q ) [l.-.5(e (6,8) 

+ £ (6,10))]} + CL. { (B -B ) * 

WC 18 6 

[1-.5E (6,8) ] } + CL (1-CL-) { (B -B ) * 

wc 1 2 10 8 

[1-.5(E (6,8) + £ (6,10))]} . (3-5) 

wc wc 
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f 4 * = f 1-CL i ] V t5[£ wc (4 ' 6)( V B 4 ) 

+ (£ wc (4 ' 6) + £ wc (4 ' 8))( VV + £ wc (2 ' 4)( V B 2 ) 

+ (£ wc (2,4) +£ wc ( 1,4))(B 2 -B 1 ) + e wc ((0,4),T 1 )*B 1 ]} 

+ CL. {B . 5 [£ (2 , 4) (B -B_) (3-6) 

1 4 WC 4 z 

+ (£ (2,4) + e (1,4))(B -B ) + e ( (0,4) ,T ) *B ] } 

WC WC Z 1 WC 1 1 

+ (1-CL 1 ) (1-CL 2 ) {(B 10 -B 8 ) [l-.5(e wc (4,8) + ^ (4 , 10) ) ] } 

F 2* = [1 " CL 1 ] {B 8"* 5[£ wc (2 ' 4)(B 4 _B 2 ) + (£ wc (2 ' 4) 

+ £ wc (2 ' 6))( VV + (£ wc (2 ' 6) + £ wc (2 ' 8)) * 

(B -B ) + E (1,2) (B -B ) + £ ( (0,2) ,T )*B ] } (3-7) 

8 6 v/c z 1 wc 1 1 

+ (l-CL ) (1-CL ) {(B. -B r ) [l-.5(E t (2,8) 

1 z 10 o WC 

+ £ wc (2 ' 10))]} + CL 1 {B 4 _ ' 5 [£ wc (2 ' 4) (B 4" B 2 ) 

+ £ (1,2) (B -B ) + £ ((0,2),T )*B ]} . 

WC Z 1 wc 1 1 

As was described by Spaeth (1975, Section III.B.5.) concerning the 

use of the composite case, Eq. (3-2) can be reduced to give expressions 
* 

for for the various reference cloud-cover cases (0,0), (1,1), (1,0) 

* * 

and (0,1). The resulting schematics in the case of F and F are 

8 4 

depicted as the unhatched area in Figs. 3(a,b,c,d) and 4(a,b,c,d), 

* * * 

respectively below. Similar graphs for F^ , F^ and F^ can be found 
on pp. 41-43 of Spaeth (1975) , and are not reproduced here. 

* 

A typical gridpoint listing of the IR net-flux computations, F 
(k = 10, 8,... 2), has been reproduced in Table II for gridpoint (1,1) 
based upon the radiative sounding of 16 April 1974 (see Table I) . The 
printout procedure involves computation of the reference net-flux 
values 

F k *(0,0), F k *(l,0), F k *(0,l), F k *(l,l) 
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0 0 
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Figure 3 C Terrestrial net flux F 8 * with (a) clear skies (case (0,0)); (b) high overcast only 
case (1,0)); (c) low overcast only (case (0,1)); (d) both high and low overcast 
(case (1,1) ) 0 Unhatched area in each case depicts F *«, 



O’O I 1 ^ ^ ^ ^ ^ N ^ ] o’O 
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Figure 4. Terrestrial net flux F * with (a) clear skies (case ),0); (b) high overcast only 
(case (1,0)); (c) low overcast only (case (0,1)); (d) both high and low overcast 
(case (1,1)). Unhatched area in each case depicts F *„ 




from each of Eqs. (3-2), (3-4), (3-5), (3-6) and (3-7) for level 

k = 10, 8,... 2 respectively. Then Eq. (3-1) with the appropriate 

weight-factors of Eqs. (2-3a,b,c,d) has been utilized to derive the 
* 

composite F (CL , CL ) values that are listed on the bottom line of 

jC -L ^ 

Table II. 



TABLE II. A sample listed of IR net-flux computations, weighting 
factors and composite values, F * (CL , CL 2 ) , as computed for grid- 
point (1,1) for 16 April 1974. K Net flux values in ly min'” 1 . 



Cloud Case 


Weight 














(CLj^ CL 2 ) 


W(CL 1 , 


CL 2 ) 


* 

F io 


* 

F 8 


* 

F 6 


* 

F 

4 


* 

F 2 


(0,0) 


.1002 




.1664 


.2288 


.2649 


.3160 


.3487 


(1,0) 


.1497 




.0943 


.1331 


.0915 


.1706 


.2338 


(0,1) 


.3007 




.0512 


.1563 


.2341 


.2871 


.3202 


(1,1) 


.4494 




.0512 


.0586 


.0607 


.1706 


.2338 


* 

F v -composite values 


.0692 


.1171 


.1379 


.2202 


.2713 


C. TROPOSPHERIC COOLING BY LAYERS 


* 

? F io 


COOLING 






For each 


mid-seasonal ■ 


day listed and at each gridpoint, 


an IR net 


flux computation in 


the format of 


Table : 


El is easily 


converted into 


four sets of 


layer cooling 


effects 


• 












* 

F810 = F 

o 


* 

- F 

10 












F68 


* 

= F 6 


* 

‘ F 8 












F46 


* 

= F 4 


* 

~ F 6 












F24 


* 

= F 2 


* 

- F 4 






(3-8) 


These layer 


cooling 


rates 


(ly min 


^ have 


then been collected 


in meri- 



dional cross-section format for each longitude under study (see Section 
V) and by season. 
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The overall tropospheric cooling rate by IR net flux is then given 
* * 

simply by - F^ Q at each gridpoint for the date of the radiational 

sounding. The tropospheric cooling rates computed are then identical 
to 

★ * 

F 2 - F = F24 + F46 + F68 + F810 . (3-9) 

★ * 

The values of - F^ so deduced are discussed on both a seasonal 

and a zonally - averaged basis in Section V. 

It will suffice to discuss here the zonally-averaged values of 

* 

F (CL^, CL^) as computed by the long-wave radiational model previously 

* 

presented in this section Q These results, listed simply as F in 
Table III, will be discussed as a function of seasons, latitudes and 

* 

CL (total opaque cloud cover given by Eq. (2-4)). The listings of F^ 

in Table III are essentially as extracted from computations in the 

* 

format of Table II followed by meridional-averaging of F across 

constant latitude lines. Finally the zonally-averaged annual values 
* 

of both F^q and CL have been computed by arithmetic-averaging over 
the four mid-seasonal results at each five-degree increment of latitude 
from 20S to 65N (cf., Eq. (5-1)). 

* 

The model-annual values of F^ presented in Table III are presented 

with those derived from Budyko (1956) , which in turn are listed in the 

final column of Table III. Corresponding values of total opaque cloud 

cover, CL, for the Budyko climatology were not available so that only 

★ 

a general comparison of the two annual F zonally-averaged distribu- 
tions is possible. 
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Annual 

16 January 16 April 16 July 16 October Model Values Budyko 
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* 

Table III depicts the zonally-distributed values of and of CL 

at the earth's surface. In a seasonal comparison of the model-computed 

* * 

values it is clearly shown that (CL) i^s ci decreasing function 

of cloud cover . There is a clear-cut tendency in each season for a 

★ 

maximum value of F^ Q to be located in the subtropics (latitudes 15N- 

* 

25N) . Also there is evidence of a high latitude (55N-60N) minimum F^ 
associated with a concentration of maximum cloud cover CL. An outstand- 
ing variation is the transition in the Southern Hemisphere latitudes 
(20S - 10S) , which has small cloud cover in local winter and compari- 
tively large cloud cover during the other three data periods. This 

cloud-cover variation corresponds in general to the ITCZ behavior in 

* 

these latitudes across the indicated seasons; so that, F^ is a maximum 

in mid-July and a relative minimum in the period January- April . 

The final entry in Table III, namely "Wt. Avg." , is the Northern 

Hemisphere mean, cosine-weighted with respect to latitude (cf., Eq. (6-1)). 

* 

The cosine-weighted F values show a minimum in Northern Hemisphere 
summer with no clear-cut differences in CL. This summer minimum is 
attributable to higher downward flux with increased summer vapor 
pressures . 

D. OUTGOING IR NET FLUX TO SPACE 

1 . Parameterization Formula for Top-of-Atmosphere IR Net Flux 

In the four previous radiation studies (A, B, C, D) , an approxi- 
mation to the IR net flux to space, designated as FF2 , was computed 

from the radiative sounding at each gridpoint. FF2 was essentially 

★ 

an extrapolation of F^ to the top of the atmosphere obtained by delet- 
ing the downward IR flux due to stratospheric water-vapor and CO^. A 
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more precise expression for the net IR flux to space was introduced here 
after Martin (1975) o The development is analogous to the quadrature 
summations 



B=0 




for the various reference-cloud combinations (0,0), (1,0), (0,1) and (1,1) 

of Fig. 5(a,b,c,d) respectively . The final quadrature-formula is then 
obtained as the weighted net-flux result as was also done in Eqs. (3-4) , 
(3-5) , (3-6) , (3-7) and is listed below. 

F T = [1~CL 1 ] {Bq-, 5 [ (^ wc (Of 1) + e wc (0'2)) (B 2 - Bl ) 

+ (£ (0,2) + £ (0,4) ) (B -B ) 

wc wc 42 

+ (£ (0,4) + e (0,6)) (B-B.) + (e (0,6) + £ (0,8)) (B -B )]} 

wc wc 6 4 wc wc 8 6 

+ CL {B-.5[(£ (0,1) + £ (0,2) ) (B-B ) (3-11) 

1 4 wc wc 2 1 

+ (e (0,2) + £ (0,4) ) (B -B ) ]} 

wc wc 4 2 

+ (1-CL ) (1-CL ) { [1- . 5 (E (0, 8) +£ (0,10))] [B-B ]} 

1 2 wc wc 10 o 

where CL^ , CL^ are given by Eq. (2-2) . Note that of Eq. (3-11) is 
representable by the unhatched areas of Figs. 5(a,b,c,d) for the various 
reference-cloud cases and that F^ has no downward IR flux corresponding 
to the level k = 0. 

2 . Comparisons of F ^ with Mid-Seasonal Satellite Climatology 

Comparison was made of computed-model values of F^ with satellite 
measurements of total long-wave flux to space (after Raschke et al., 1973) 
for the Nimbus III mid-seasonal periods most nearly comparable to that 
of the FNWC data. F^ for each mid-seasonal date was computed over the 
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Figure 5. Outgoing terrestrial net flux to space with (a) clear skies [case (0,0)]; (b) high ove 

cast only [case (1,0)]; (c) low overcast only (0,1); (d) both high and low overcast [case 

(1,1)]. Unhatched area in each case depicts F . 




four meridians considered. These F t values were then averaged across 
the four meridians to get a mean zonal distribution of the type shown 
in Table IV. 

Table IV shows the zonally-averaged model values compared with 
those extracted from Raschke. Raschke's results were obtained by 
averaging across the same oceanic meridians as those used in this study. 
Again, in the bottom line of each column in Table IV is listed the 
cosine-weighted mean of each set of column values for the Northern 
Hemisphere only. 

The zonally-averaged values computed by the F^-model are very 
close to those reported by Raschke, especially between 0-65N and in 
the Northern Hemispheric means . The limitations of the comparisons 
made here are obvious, when it is recalled that between latitudes 20S- 
5S and between 60N-65N there are fewer than four meridional lines 
available for computing the listed zonal values in Table IV. For all 
other zonally-averaged values, four meridional lines were used in the 
averaging. 

The close comparison between model-values of F^ from (3-11) and those 
essentially derived from satellite climatology tend to support the 
cloud parameterization, Eq. (2-2) , insofar as IR net flux is concerned. 
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IV. SOLAR RADIATION 



A. COMPOSITION OF SOLAR INSOLATION 

At the top of the atmosphere (k=0) this study assumed a solar con- 
stant of 2.00 ly min 1 (Joseph, 1971). Furthermore, this constant was 
assumed subject to a four percent attenuation above the tropopause due 
to ozone and oxygen. Thus the effective solar constant at level k=2 
in this study is 1.92 ly min 1 . 

To compute the effective solar insolation at the tropopause the 
following formula was used 

F (2) = S [— — ] ” 2 Cos Z (4-1) 

m 

where S is the effective solar constant at level k=2 and 

r/r = ratio of the actual earth-sun distance to the mean 
m 

earth-sun distance, a function of the Julian date 
Cos Z= cosine of the zenith angle, a function of the 
Julian date determined by 

Cos Z = Sin <j) Sin 6 + Cos <j) Cos 6 Cos h . (4-2) 

The symbols on the right side of Eq. (4-2) are defined as follows: 

(J) = latitude 
6 = solar declination 
h = hour angle 
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The Smithsonian Meteorological Tables (List, 1958) list the ratio of 



r/r and the solar declination, 6, for the mid-seasonal dates appli- 
m 

cable to this study and reproduced in Table V. 



TABLE V. Values of the ratio of the earth-sun radius vector, r/r , 
and of the solar declination angle, 6, used in this study for year 



1973 - 1974. 






Date 


r/r 

m 


6 


16 January 


0.98372 


21 . 07917°S 


16 April 


1.00333 


8 . 48333°N 


16 July 


1.01644 


21 . 50000°N 


16 October 


0.99717 


8. 22500°S 



The value of sin (p was calculated using one of two different formulas, 
depending on the data-line used for the computations, in terms of the 
FNWC map coordinates (I,J) as in Eq. (4-3a,b) . Conversely for these 
lines one may solve for I in terms of sin (p as in (4-3c,d) : 



Lines 1,3,4 



Line 2 



Lines 1,3,4 



Sin (p 



Sin (p 



973.752 - 2 (32-1) 2 
973.752 + 2 (32-1) 2 
973.752 - (32-1) 2 
973.752 + (32-1) 2 



I = 32 - 22.065 [ 



Cos (p 
1 + Sin (p 



(4-3a) 



(4-3b) 



(4-3c) 



Line 2 




31.205 

,25 

,25 

,38 



Cos (p 
1 + Sin <J) 

for Lines 1,2 
for Line 3 
for Line 4 



(4- 3d) 
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Here I is the abscissa distance on the FNWC grid (Fig. 1) and varies 
by line as described in Section II. The soundings for lines 1, 2, and 
3 were all taken at OOOOGMT with the solar noon existing at the 180th 
meridian; therefore the hour angles for these three lines were 55°, 

10° , and 35° , respectively. For line 4, the soundings were taken 12 
hours earlier with solar noon at the Greenwich meridian, giving an 
hour angle for line 4 of 35°. 

A very simple partition of solar insolation was utilized in this 
study after Joseph (1971) . It consisted of dividing the insolation 
F(2) into two parts at level k=2. One part was considered to include 
all wavelengths X > . 9 ym where absorption by water vapor and carbon 
dioxide bands are the most prevalent attenuation processes in clear air. 
This part of the solar spectrum was termed the F (A) energy and con- 
sidered subject to water-vapor absorption but not to Rayleigh scatter- 
ing. For shorter wavelengths, X < .9 ym, absorption of the solar in- 
solation energy by water vapor was considered negligible. This part of 
the solar insolation was denoted F(S) suggestive of the fact that it 
was subject only to Rayleigh scattering attenuation in clear air. The 
two partitions are formulated after Joseph (1971) as follows: 



F (A) 


= .349 F (2) 


(4-4) 


F (S) 


= .651 F (2) . 


(4-5) 



In this study, the introduction of two cloud decks produced cloud-reflec- 
tivity effects upon both the F (A) and F(S) solar energy insolations. 
However, in the clear areas around any gridpoint only the absorption- 
attenuation applies to the F (A) insolation, while only Rayleigh scatter- 
ing-attenuation applies to the F(S) insolation. 
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B. DISPOSITION OF F(S) INSOLATION 



In the disposition of the F(S) insolation, Joseph (1971) determined 
that Rayleigh scattering reflectance to space by clear skies (after 
Coulson, 1959) could be effectively approximated by least squares in 
the following form 

a (R) = .085 + .25074 [log (^-) Sec z) ] (4-6) 

o 

where P = 1013.25 mb. In Eq. (4-6), tt/P = 1 in view of the fact that 
o o 

mean sea level pressure 7T is close to 1000 mb. Also 

Sec z = (Cos z) 

with Cos z given by Eq. (4-2). 

The surface albedo a(G) is another reflective parameter utilized in 
this study. Over oceanic areas the following formula for a(G) after 
Gates et al (1971) , was utilized: 

a(G) = max {.06, .06 + .54 (.7 - Cos z) } . (4-7) 

As described in Section III, four combinations of reference-cloud 
cases are possible with a two- layer cloud model. The disposition of 
F(S) under each of these cases will be discussed in the remainder of 
this subsection. 

1. Clear Sky Case 

In the clear sky (0,0) case the F(S) insolation was subjected 
to both Rayleigh scattering reflectance a(R) and the surface reflectance 
a(G) . Considering the likelihood of a succession of multiple reflections 
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between earth and atmosphere, the F(S) insolation actually penetrating 



the earth’s surface after scattering is given by 



IS10(0,0) = F(S) [l-a(R)] [l+a(R)a(G) + ... (a(R)a(G) ) 



n 



+ . ..] * (l-a(G) ) 



(4-8a) 



that is, by 



isio(o,0) = F (s ) [l-a(R) ] [i-a(G)]/[i-a(R)a(G) ] 



(4-8b) 



2 . Cloudy-Sky Cases 

In the three cases in which clouds were present, F(S) insolation 
absorbed by the ground at each gridpoint was computed using the follow- 
ing equation (after Arakawa, 1972) : 



As indicated by the notation (1,1), denoted CL^ = CL^ = 1.0, Eq. (4-9) 
is the formula used in calculating F(S) insolation absorbed by the 
ground in the case where overcast clouds are present at both levels of 
Fig. 2. Also in Eq. (4-9), initial values of cloud-reflectance were 
chosen, namely R(l) = .54 for the mid-level clouds between k=4 and 6, 
and R (2 ) = .66 for the low-level clouds between k=8 and 9. Both cloud- 
reflectance values are as suggested by C. D. Rodgers (1967) . 

For all the other cloud cases, the following changes were 
applied to Eq. (4-9). In the (1,0) case (CL^ = 1.0, CL 2 = 0.0), the 



IS10 (1 , 1) = F (S ) (l-R(l)) (1-R(2) ) (l-a(G)) 



*{l-[R(l)R(2) + R(2)a(G) + R(l)a (G) 
- 2R(l)R(2)a(G)] } _1 . 



(4-9) 
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desired earth-absorbed insolation is obtained by setting R(2) =0.0 
in (4-9) , from which it follows that 

IS10 ( 1 , 0) = F(S) (l-R(l) ) (l-a(G) )/[l-R(l)a(G) ] . (4-10) 

In the case (0,1), one sets R(l) = 0.0 in (4-9) so that (4-9) simpli- 
fies to 

IS10(0,1) = F(S) (1-R (2 ) ) (l-a(G) )/[l-R(2)a(G) ] . (4-11) 

Note that with a cloud overcast present, the Rayleigh clear-sky scatter- 
ing a(R) does not appear in Eqs. (4-9), (4-10) or (4-11), but is included 

empirically in the cloud reflectances R(l) and R(2) . 

3 . Composite F(S) Insolation at Earth 

Equations (4-8) , (4-9) , (4-10) and (4-11) were utilized in the 

computation of the cloud-weighted F(S) insolation penetrating the 
earth 1 s surface considering the areal-weights of the cloud combinations 
denoted by (0,0), (1,1), (1,0) and (0,1) about a gridpoint. The re- 

sultant F (S ) insolation penetrating the earth’s surface denoted by 
IS10 is therefore expressible as 

IS10(CL 1 , CL 2 ) = IS10(0,0) W(0,0) 

+ IS10(1,1) W(l,l) 

+ IS10(1,0) W(1,0) 

+IS10(0,1) W(0,1) . (4-12) 

Here the weighting factors W(0,0), W(l,l), W(1,0) and W(0,1) are com- 
puted in Eqs. (2-3a,b,c,d) respectively. Note finally that the part 
of F (S ) insolation reflected to space is found by subtracting IS10 
(CI^, CL 2 ) from F(S) . 
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Table VI lists the F (S ) -disposition resulting from a particular radia- 
tive sounding at gridpoint (1,1) on 16 April 1974. The individual 
computations of IS10 as they apply for the possible overcast-clear 
layer cases are made under the heading "ISIO." The difference 

REFS = F(S) - ISIO (4-13) 

in each case represents F (S) -insolation reflected to space while 

STEAN ■ -iSk (4 - 14) 

has been computed as that portion of the F (S) -insolation incident at 
the sea surface just prior to transmission by the surface. Note that 
no absorption in air has been included in the computations of Table VI, 
and that the only absorption permitted is that implicit in ISIO. 

Finally at the bottom of each column, e.g., ISIO, the composite value 
has been computed by means of the weighting scheme of Eq. (4-12) . 

TABLE VI. A sample listing of values of F(S) insolation (ly min 
computed at gridpoint (1,1) for 16 April 1974 using Eqs. (4-6,..., 

4-14) . 

Cloud-case Weight 



(Cl , cl 2 ) 


W(CL lf CL 2 ) 


ISIO 


REFS 


STRAN 


(0,0) 


.0993 


.4305 


d751 


.5216 


(1,0) 


.1484 


.2539 


.3517 


.3076 


(0,1) 


.3016 


.1921 


.4135 


.2327 


(1,1) 


.4507 


.1340 


.4716 


.1696 


Composite-F (S) 


values 


.2014 


.4041 


.2441 
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C. DISPOSITION OF F (A) INSOLATION 



The fractional portion of the solar insolation subject to absorption 
by atmospheric water-vapor and carbon dioxide is covered in the follow- 
ing subsections. 

1. Clear-sky Case (0,0) 

The Manabe-Moller absorptivity function provided the necessary 
absorptivity values for the key layers in this case. The form of this 
absorptivity function is 

303 

a(2,k) = . 271 [U (2 , k) Sec z] * . (4-15) 

Here a is the absorptivity applied to the pressure-scaled water vapor 
mass between levels 2 and k (Fig. 2) along the zenith slant-path angle 
z. The resultant absorbed insolational energy in the particular layer 
(2,4) is then given by the Manabe-Moller relation 

303 

A24 = 0.271F(A) [U(2,4) Sec z] # . (4-16) 

In the same manner A26, A28 and A210 are found. Then the absorbed in- 
solation in the layers (4,6), (6,8) and (8,10) are computed by 



A46 = A26 - A24 


(4-17) 


A68 = A28 - A26 


(4-18) 


A810= A210- A28 . 


(4-19) 



Water-vapor mass above level 2 was assumed negligible in the F (A) dis- 
position of the solar insolation. 

By subtracting A210 from F(A), the direct transmission of F (A) 
insolation impinging at the earth's surface was determined. The trans- 
mission of F (A) insolation is then further reduced by the transmissivity 
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(l-a(G)) after surface-ref lectance, which leads to the earth-absorbed 
insolation 

IA10 (0,0) = F(A) {1-.271 [U(2, 10)Sec z] * 3 ° 3 } (l-a(G) ) . (4-20) 

The transmitted energy impinging upon the earth just prior to absorp- 
tion is 



TRANA (0,0) = IA10(0,0)/[l-a(G) ] . (4-21) 

In the remainder of this subsection, which discusses the cloudy 
layer cases, representative cloud reflectivities and cloud absorptivi- 
ties were initially adopted, after C. D. Rodgers (1967) , for the two 
possible cloud layers. These initial cloud reflectivities were RA(1) = 
.46 and RA(2) = .50 while the cloud absorptivities were A(l) = .20 and 
A (2) = .30. Note that the reflectivities for the F (A) wavelengths are 
somewhat smaller than those adopted for the F(S) wavelengths. The pro- 
cedure of considering the cloud conditions to be overcast whenever they 
appear and then applying the appropriate weighting factors in the com- 
posite summation will again be followed as in Sec. IV. B. 

To simplify the discussion for the cloud-covered cases, the (1,1) 
case will be presented first as it contains representative type equa- 
tions for the remaining two cases, (1,0) and (0,1). 

2 . Overcast in Both High- and Low-cloud Layers 

The following set of formulas illustrate the model disposition 
of incoming solar insolation (F(A)) from level k=2 to the earth* s sur- 
face and permits determination of the amount of insolation absorbed by 
the atmospheric layers and by the earth* s surface. The dashed separa- 
tion lines are introduced to subdivide the absorption and reflection 
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physics of the model into subsections which permit the analysis to 
proceed more or less within successive 200 mb layers 0 The equations 
relate to the parameters in Fig* 6 where the insolations , A24 , A46, 

A68, A89 and A910, etc., represent the contributions to the insola- 
tion absorbed in the layers involved. Symbols F2, F4, F6, F8 and F9, 
etc., depict the streams of insolation passing through the indicated 
level. A vertical arrow implies the direction of insolation passage, 
i.e., 1 denotes downward insolation, f upward-reflected insolation, 
and downward-reflected insolation. Terms involving the symbols 
"TD", as expressed by the functions of Eq. (4-22), indicate the Manabe- 
Moller transmissivities for diffuse insolation in the layer beneath an 
existing cloud. In the latter situation, the term Sec Z in Eq. (4-15) 
is effectively replaced by the mean slant-path, Sec Z — 5/3 (Katayama, 
1966) , e .g . , 

TD68 = 1- . 271 [U (6,8) 5/3] * 3 ° 3 

303 

TD910 = 1-.271 [U (9,10) 5/3] * . (4-22) 

In accordance with the above definitions and as described by Fig. 6, 
the formulas for the (1,1) case are listed: 

F4+ = F(A) (1-. 271 [U (2,4) Sec z]* 303 ) 

F4+ = F44' (RA (1) ) 

303 

F2+ = F4+ (1- . 271 [U (2 , 4) Sec z] ) 

A24 = F (A) - F44- + F4+ - F2+ (4-23a) 



48 



A46 = F4+(A(1)) 



(4-23b) 



F64 = 


F41 - F4t - A46 


F84 = 


F64 (TD68) 


F8t = 


F84 (RA (2) ) 


F6t = 


F8t (TD68) 


F644= 


F6t (RA (1) ) 


F844= 


F644 (TD68) 


A68 = 


F6l - F84 + F8t - F6t + Fell - F844 (4-23c) 



A89 = 


(F84 + F844) (A ( 2) ) 


F94- = 


F84 - F8t + F844 - A89 


F104= 


F94 (TD910) 


F10t= 


FI 04 (a(G) ) 


F9t = 


FlOt (TD910) 


F9i4= 


F9t (RA (2) ) 



F10H= F9H (TD910) 

A910 = F9I - F104 + FlOt - F9+ + F944 - F1044 

A810 = A89 + A91C (4-23d) 

Note that the effect of multiple reflections between clouds or 
between the earth's surface and the lower cloud has been incorporated 
to include only the effect of two reflections, with the lowermost re- 
flecting surface absorbing the remaining impinging insolation. Computa 
tions indicated that the insolation remaining after two reflections was 
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Figure 6. Schematic representation of F(A) insolation disposition in 
the case of two overcast layers. 
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too small to warrant the consideration of further reflections. Also 
insolation reflected upward from a lower interface (cloud or ground) 
to the base of an upper cloud deck has not been subjected to upper 
cloud absorption. This tends to reduce very slightly the secondary 
cloud-absorption. 

From Eq. (4-23d) , the impinging F (A) insolation at the earth’s 
surface can be expressed 



The F (A) insolation which is actually absorbed by the earth's surface 
(see Fig. 6) may be written as 



3 . Disposition of F(A) Insolation with an Upper Overcast Only 

With a single cloud layer present only at the upper level (Fig. 
7) the equations depicting the model-disposition of incoming insolation 
becomes a simplified subset of the previous case: 



TRANA (1,1) = FlOl + F10+ + 



(4-24) 



IA10 (1,1) = F10l(l-a(G)) + F10++ 



(4-25) 



F4+ = F (A) (1- . 271 [U ( 2 , 4) Sec z] 



.303 



F4t = F4+ (RA ( 1) ) 

F2+ = F4t(l-.271[U<2,4) Sec z]* 303 



A24 = F (A) - F4+ + F4i - F2t 



(4-26a) 



A46 = F44- (Al) ) 



(4-26b) 
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Figure 7. 



Schematic representation at F(A) insolation disposition 
with an upper overcast layer only. 
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F64- 


F44- - F4t - A46 


F84- 


F64- (TD68) 


F104- = 


F64- (TD610) 


F lOt = 


F104- (a(G) ) 


F8+ 


FlOt (TD810) 


F6f 


F10I (TD610) 


F64-+ = 


F6t (RA ( 1) ) 


F84-4- = 


F64-4- (TD68) 


FlOii = 


F84-4- (TD810) 


A68 


F64- - F8+- + F8t - F6+ + F6ll - F844 (4-26c) 



A810 = 


F84- - F10+ + FlOf - F8+ + F844 - F1044 (4-26d) 



The variables used above are defined in a similar manner to those in the 
(1,1) case* The impinging insolation at the earth’s surface may be for- 
mulated as follows 

TRANA (1,0) = F10+ + FlOK (4-27) 

while the F(A) insolation absorbed by the earth in this case is given by 
IA10 (1,0) = Fl0i(l-a(G)) + F1044 . (4-28) 
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4. Disposition of F (A) Insolation with a Low Overcast Only 



With an overcast lower cloud layer the model-disposition symbols 
are as depicted in Fig. 8 and are physically related as follows: 



F44- 


= F (A) (1- . 271 [U (2 ,4) 


_ , • 303 

Sec z] ) 




F64- 


= F (A) (1- . 271 [U (2 , 6) 


_ , - 303 

Sec z] ) 




F84- 


= F (A) (1- . 271 [U ( 2 , 8) 


o , .303, 

Sec z] ) 




F8t 


= F8HRA(2) ) 






F6t 


= F8t(l-.271 [U (6 , 8) 


, .303, 
Sec z] ) 




F4 + 


= F8+ (1- . 271 [U (4 ,8) 


„ ,-303, 

Sec z] ) 




F2 + 


= F8t (1- . 271 [U(2,8) 


„ ,.303, 

Sec z] ) 




A24 


= F (A) - F44- + F4 + - 


F2 + 


(4-29a) 


A46 


= F44- - F6+ + F6+ - 


F4 + 


(4-29b) 


A68 


= F64- - F84- + F8I - 


F6t 


(4-29c) 


A89 


= F84- (A ( 2) ) 






F94- 


= F84- - F8+ - A89 






F104- 


= F94- (TD910) 






F10+ 


= Fl(H (a (G) ) 






F9i 


= fio-Mtosio) 






F94-4- 


= F9+ (RA(2) ) 






FlO'1'4- 


= F94-4- (TD910) 






A910 


= F94- - F104- + FlOt 


- F9+ + F94-4- - F1044 




A810 


= A89 + A910 




(4-29d) 
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Figure 8. Schematic representation of F(A) insolation disposition 
with a lower overcast layer only. 
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Again the variables are defined as stated before in the (1,1) 



case. Likewise the incident flux is defined at the earth* s surface as 



In Eqs. (4-25) , (4-28) , and (4-31) , the quantity FlOil is small enough 
in each case so that no further reflections from the earth were con- 
sidered. 

5. Composite F(A) Layer-Absorptions and Surface-Absorption 
Insolation 

As has been previously discussed, the standard grid-area weight- 
ing scheme of this study was applied to obtain composite values of the 
absorbed F (A) -insolation in key layers and also within the earth * s sur- 
face. The weighting factors applied to the corresponding overcast- 
combination absorption quantities provided the following composite 
results : 



TRANA (0,1) = FlOl + FlOll 



(4-30) 



while that portion which is absorbed by the earth (Fig. 8) is 



IA10 (0,1) = Fl0i(l-a(G)) + F1011 



(4-31) 



A24(CL 1 ,CL 2 ) = A24 (0,0) W(0,0) + A24(l,0) W(1,0) 



+ A24 (0,1) W(0,1) + A24 (1,1) W(l,l) 



(4-32) 



A46(CL 1 ,CL 2 ) = A46 (0,0) W(0,0) + A46(l,0) W(1,0) 



+ A46 (0,1) W(0,1) + A46 (1,1) W(l,l) 



(4-33) 



A46(CL l ,CL 2 ) = A68 (0 , 0) W(0,0) + A68(l,0) W(1,0) 



+ A68 (0,1) W(0,1) + A68 (1,1) W(l,l) 



(4-34) 
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(4-35 ) 



A810(CL 1/ CL 2 ) = A810 (0, 0) W(0,0) + A810(l,0) W(1,0) 

+ A810 (0,1) W(0,1) + A810) 1,1) W(l,l) 

IA10 (CL 1 , CL 0 ) = IA10 (0,0) W(0,0) + IA10(1,0) W(1,0) 

+ IA10(0,1) W (o , 1) + IA10 (1,1) W(l,l) • (4-36) 

The weighting factors W (0,0) , . • . , (Wl, 1) were listed in Eqs . (2-3a ,b, c ,d) , 

and A24 (0,0) , A46(0 / 0) / A68(0,0), A810(0,0) and IA10(0,0) are given in 
each clear sky case (0,0) about each gridpoint by Eqs. (4-15) , (4-16) , (4-17) 
(4-18) , (4-19) and (4-20) respectively. 

The results for the absorption in layers (2,4) , (4,6) , (6,8) , 

(8,10) and at the surface, level k = 10, are shown in the following table: 

TABLE VII. A sample listed of F (A) insolation values (ly min 1 ) 



computed 


at gridpoint 


(1,1) for 


the 16 


April 


case. 








Cloud-Case Weight 

CL , CL 2 W(CL lf CL 2 ) 


A24 


A46 


A68 


A810 


IA10 


REFA 


TRANA 


(0,0) 


.0993 


.0454 


.0356 


.0386 


.0418 


.1348 


.0285 


.1634 


(1,0) 


.1484 


.0633 


.0559 


.0280 


.0156 


.0482 


.1137 


.0581 


(0,1) 


.3016 


.0464 


.0381 


.0728 


.0771 


.0239 


.0664 


.0286 


(1,1) 


.4507 


.0633 


.0559 


.0432 


.0285 


.0104 


.1233 


.0126 


Composite 


-F (A) values 


.0564 


.0485 


.0494 


.0425 


.0324 


.0953 


.0391 



In the computational scheme indicated by the entries of Table 
VII, the reflected F (A) insolation to space has been depicted by the 
symbol REFA, and its values follow from 

REFA = F (A) - A24 - A46 - A68 - A810 - IA10 (4-37) 
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whereas the TRANA dispositions are given by Eqs . (4-21), (4-24), (4-27) 



and (4-30) or by its weighted-mean value in the case of TRANA-composite . 

6, Absorptivity (ABA) by Layers 

Here the (fractional) absorptivity as well as the actual insola- 
tion values absorbed in the layers are considered. In the computation 
of absorptivity, which is fractional absorption, the total undepleted 
insolation at the top (k=0) is used as a base. The following equation 
was utilized in this calculation: 

-2 

FADJ = 2.00 (r/r ) Cos z . (4-38) 

m 

The absorptivity of the troposphere ABA was computed from the ratio of 
the insolation absorbed in the troposphere to the insolation incident 
at the top of the atmosphere rather than at k=2 : 



ABA 



A24 + A46 + A68 + A810 
FADJ 



(4-39) 



D. ALBEDO (ALB) OF THE EARTH -TROPOSPHERE SYSTEM 

In considering the planetary albedo, the reflected insolations of 
the earth-troposphere systems in both the F(A) and F(S) insolational 
regions must be recalled by the program. Thus REF is computed at each 
gridpoint as the sum of the reflected insolation energy in F(A), denoted 
REFA in Eq. (4-37) and the reflected part of F(S) denoted REFS in Eq. 
(4-13) : 



REF = REFS + REFA (4-40) 

Finally the planetary albedo is related to FADJ through 
REF 



ALB 



'FADJ # 



(4-41) 



E. COMPOSITE ABSORPTIVITY (ABG) BY THE EARTH- SURFACE ; 

COMPOSITE ATMOSPHERIC TRANSMISSIVITY (ATRAN) 

1. Absorptivity (ABG) of Earth 

By summing the weighted values of F(S) and F(A) portions of 
the incoming insolation entering the earth, the total insolation ab- 
sorbed at the earth's surface was computed. This quantity when 
divided by the extraterrestrial insolation gave the fractional absorp- 
tivity (ABG) of the earth's surface. The equation for ABG was 



ABG = 



IA10 + IS10 
FADJ 



(4-42) 



where IA10, IS10, and FADJ were defined previously by Eqs. (4-36), 
(4-12) and (4-38) respectively. 

2 . Transmissivity (ATRAN) of the Troposphere 

Also computed was the total insolational energy TRAN, incident 
at the earth's surface just before absorption by the surface. This 
calculation is given by 



TRAN = TRANA + STRAN . 



(4-43) 



Here STRAN = [IS10/ (1-a (G) ] was previously defined in Eq. (4-14) and 
a(G) was given in Eq. (4-7) . TRANA has also been defined as the weighted 
value of TRANA (0,0), TRANA (1,1), TRANA (1,0) and TRANA (0,1) given by 
Eqs. (4-21), (4-24), (4-27) and (4-30). Note also in justification of 

STRAN that the four cases for IS10 of (4-8), (4-9), (4-10), and (4-11) 

each have the common factor (l-a(G)) in the numerator and therefore each 
transmitted F(S) insolation component available at the earth just before 
absorption needs only be divided by (l-a(G)). TRAN may thus be viewed 
as the total insolational energy incident at a pyrheliometer located at 
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earth. The (fractional) transmissivity of the troposphere (ATRAN) is 



then computed from 



ATRAN E TRAN/FAD J . (4-44) 

Note finally that the major dispositions of the total insolation at 
the indicated map times have now been identified by the fractional 
values, ALB, ABA or ABG, and ATRAN, representing the reflectivity 
(albedo) , absorptivity of air or earth, or atmospheric transmissivity 
as the case may be. 

F. ALBEDO TUNING BY COMPARISONS WITH SATELLITE CLIMATOLOGY 
1 . General Remarks Concerning a Need for Tuning Albedos 

In the four previous radiational studies (A, B, C, D) compari- 
sons were made between the solar-insolation albedo-model (ALB MOD) com- 
putations and the satellite-climatology albedo (ALBRAS) of Raschke 
et al. (1973). The results of these comparisons indicated excessively 
high values of ALBMOD especially in the tropical and subtropical oceanic 
areas. The question was raised by the previous investigators as to 
whether vertically-structured convective cloud elements in the tropics 
would have as high a reflective capability as attributed to the large- 
scale cloud masses existing primarily in horizontal layers as indicated 
by the parameterization of Eq. (2-2) . It should be recalled that the 
initial cloud-reflectances used were as suggested by C. D. Rodgers 
(1967) and are listed in Secs. IV. B. 2. and IV.C.l. 

For the purpose of tuning the model solar dispositions, and 
in particular ALBMOD with ALBRAS, the initial values of Rodgers* cloud- 
reflectances were simply adjusted by a multiplicative factor (f) which 
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turns out by the least-squares fitting to be smaller than unity. This 
device allows a greater fraction of solar radiation to penetrate down- 
ward through the cloud-types in both the tropical and subtropical areas 
and provides better agreement between ALBMOD and ALBRAS. 

2 . Method of Tuning ALBMOD to ALBRAS 

The model- tuning process involved several steps. First, the 
new variable f was defined and allowed to range in value from 1.2 to 
0.25 by 0.05 intervals. The f-values were then multiplied by the 
initially used cloud-reflectance values (after Rodgers, 1967) giving 
twenty sets of new cloud-reflectivity values. For each set of new 
values, ALBMOD and subsequently the difference A (ALB) between ALBMOD 
and ALBRAS ; i . e . , 



A (ALB) = ALBMOD - ALBRAS (4-45) 

was recomputed at each gridpoint for the four midseasonal dates. 

Because cloud-types differ in structure over tropical and extra 
tropical oceanic areas the resulting A(ALB) data, from (4-45), were 
divided into two populations denoted TR and HL. The population TR in- 
cluded the set of gridpoints located south of and including 25N, while 
HL included the set of gridpoints located north of 25N over all four 
meridians . 

Only those gridpoints where the value of at least one of CL^ 
and CL^ was greater than 0.1 were included in those populations. The 
gridpoints so excluded were considered climatologically unrepresenta- 
tive for the purpose of cloud- tuning the albedo. The A (ALB) data was 
additionally divided by season giving a total of eight separate cases 
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further defined in sequence as WIN-TR, WIN-HL, SPR-TR, SPR-HL, SUM-TR, 



SUM-HL, AUT-TR and AUT-HL. 

The Bimedical set of programs (Dixon, 1973) was utilized to 
compute the root-mean-square (RMS) of A (ALB) values corresponding to 
a wide range of f-values for each of the seasons and regions. Values 
of f leading to the minimum RMS of A (ALB) were isolated and appear in 
Table VIII, along with the RMS difference of A (ALB). 



TABLE VIII. Values of f and corresponding minimum RMS 
values of A (ALB). 





f 


Minimum RMS 
of A (ALB) 


WIN-TR 


.40 


.069 


WIN-HL 


• 

00 

o 


.073 


SPR-TR 


.35 


.067 


SPR-HL 


.65 


.082 


SUM-TR 


.60 


.079 


SUM-HL 


.70 


.109 


AUT-TR 


.40 


.052 


AUT-HL 


.60 


.073 



As can be seen from Table VIII, the listed values of f in 
winter indicates the largest difference between the tropical and extra- 
tropical cloud- types. The evenly layered cloud formations of the 
northern latitudes in mid-January seems more closely to approximate 
those considered by C. D. Rodgers (1967). On the other hand, the 
smallest difference in f-values occurs between SUM-TR and SUM-HL 
indicating a less significant difference in cloud types over the 
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complete latitude range during mid-July. As would be expected the 
Spring and Autumn f-ranges fall in between the Winter and Summer cases. 

A full listing of zonally-averaged ALBMOD values, after tuning at all 
gridpoints, and of zonally-averaged ALBRAS values, generated from 
identical gridpoints, can be found in Table IX. 

Upon examination of Table IX it is apparent that the tuned 
model-albedo values compare closely with the satellite albedos derived 
from Raschke. Table IX includes the annual averages and the Northern 
Hemisphere cosine-weighted means as well as the mid-seasonal averages. 

Further inspection of the weighted means shows that ALBMOD is 
somewhat smaller than ALBRAS for each of the four seasons and annual 
average. However most of this underestimate seems to occur in the 
HL-population where the result of minimizing A (ALB) led to values of f 
substantially less than unity (Table VIII) . This result was unexpected 
and is probably due to a mismatch between the model clouds (Eq. (2-2)) 
and those averaged in the satellite climatology particularly in the 
higher latitudes. 

The conclusion is that reflectance-tuning to satellite albedo 
appears valid in the tropics and subtropics where cloud-formations 
on a given day tend to be persistent. However tuning by this method in 
high latitudes is less conclusive because interdiurnal cloud variability 
is much greater. 

In any case, all further computations of solar dispositions in 
this study will utilize the appropriate cloud-reflectance tuning factors 
as listed in Table VIII, depending upon season and latitude. 
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WINTER SPRING SUMMER AUTUMN ANNUAL 
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V. MERIDIONAL CROSS-SECTIONAL DEPICTION 
OF THE RADIATIVE-BALANCE COMPUTATIONS 



A. GENERAL 

The general design of this section is to utilize all of the compu- 
tational concepts discussed in Sections III and IV in the computations 
for a single time-step in the radiative heating package developed for 
use in the FNWC prediction model. After testing ALBMOD with the corres- 
ponding values of Raschke et al. (1973) , it was decided that only the 
computations by the appropriate cloud-reflectance tuning factors (Table 
VIII) would be displayed in the meridional cross-sections (Figs. 10... 
17) and in the comparisons of the four-layer versus the two-layer cool- 
ing rates summarized in Table X. 

The radiative calculations were performed at each gridpoint of the 
four meridians for each of the four mid-seasonal days considered; how- 
ever, only the winter and summer cases, as illustrating the two extreme 
situations, are presented in this section. 

B. GEOGRAPHICAL REPRESENTATION OF THE RADIATIVE -BALANCE DISTRIBUTION 

The FNWC gridpoint processed analyses for 0000GMT were used at the 

three Pacific cross-sections, while that for 1200GMT were used for the 
single Atlantic meridian. This was done so that the set of gridpoints 
was considered to be subject to the actual radiative- transfer calcula- 
tions involved for these specific times. Figure 9 depicts in symbolic 
language the key to the computational entries in Figs. 10, 11,... 16, 

17. This symbolic list presents the computations made at each radia- 
tive sounding gridpoint ( I, J) having data in the form of Table I. The 
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computations proceed from the tropopause (approximately level k=2) to 
the ocean surface. For purposes of climatological data comparison, 

Figs, 10, 11,. ..16, 17 were developed by interpolating gridpoint results 
to integral multiples of 5-degree latitudinal increments. The inter- 
polation routine to this gridpoint spacing made use of the Lagrangian 
cubic interpolation scheme (after Spaeth, 1975) . 



Q(D = 



(1-2) (1-3) (1-4) 
(1-2) (1-3) (1-4) 



O (I-D d-3) (1-4) 
y 2 (2-1) (2-3) (2-4) 



(1-1) (1-2) (1-4) 
y 3 (3-1) (3-2) (3-4) 



(1-1) (1-2) (1-3) 
~4 (4-1) (4-2) (4-3) 



(5-1) 



Finally for ease in reconciling the magnitudes of all radiative-trans- 
fer rates, the time-dependent solar disposition rates have been 
averaged to 24-hourly rates . 



C. EXPLANATION OF SYMBOLIC TERMS 
1 . Cross-Section at Level k=2 

The discussion of all insolation parameters discussed previous- 
ly in Section IV dealt with the specific time of day that corresponded 
to the hour angle h for the instantaneous time t under consideration. 
The incident solar insolation dealt with is then 



F (2) = S(— ) 2 Cos z . 
r 

m 



(5-2) 



In order to avoid reference to specific map times t, the instantaneous 
solar hour-angles were h = 35, 10, 55, and 35 degrees, respectively 
for cross-sections 1,2,3 and 4 as depicted in Figs. 10, 11, 12 and 13 
for winter, and in Figs. 14, 15, 16, and 17 for summer. 
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QAVE represents the 24-hour average of F(2) and appears as the 
first input symbol in Fig. 9. Its value is considered to be more re- 
presentative climatologically for the data day under consideration than 
F (2) . 

QAVE is derived by the formula 

« RVE ■ F<2) §sri <5 - 3) 



where 

Cos z = [H Sin(f)Sin6 + Cos4>Cos6SinH] /TT (5-4) 

H = ArcCos [-Tan<{>Tan<5] (5-5) 

Here 6 is the appropriate solar declination angle as listed in Table V, 
and H is the appropriate hour angle at local sunset at latitude (j). The 
value of H also depends upon which mid-seasonal date is being considered. 
Cos z in Eq. (5-4) is equal to the 24-hour average cosine of the zenith 
angle, Eq. (4-2) . The 24-hour time averaging period for QAVE gives heat- 
ing results consistent in magnitude with the terrestrial flux divergences, 
which change only slightly with the time of day. The conversion to ex- 
pected daily averaged solar disposition quantities is compatible with 
the determination of a hemispheric radiative balance for the given mid- 
seasonal dates (Sec. IV. A.). 

Other parameters needed for level k=2 are 

pr-) c 7 

QREF = REF (t) ( ) (5-6) 

Cos 2 

where 

REF (t) = F (2) - A24 - A46 - A68 - A810 

- (IA10 + IS10) (5-7) 
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REF ( t) is the instantaneous solar reflected insolation at a gridpoint 
and QREF is its 24-hour average, assuming that the instantaneous 
planetary albedo remains constant for the 24-hour period. This assump- 
tion requires that the cloud amounts computed at the indicated synoptic 
times are representative of the entire day. 

The same principle will be used with regard to all other solar 
parameters in the conversion from time-dependent values at solar time t 
to 24-hour averaged values. Superior bars ( ) are not used in the sym- 

bolism for the averaged values shown in the cross-sections key. Fig. 9, 
but are implied by the use of QAVE, etc. in the solar-disposition terms. 
The 24-hour average tropopause balance, BALT, is computed from 

BALT = QAVE - (QREF + F 2 *) (5-8) 

for the level k=2 at the indicated latitude. Net terrestrial fluxes, 

* 

such as F^ , were considered to be constant throughout the 24-hour period, 
a valid assumption if the cloud cover remains quasi-constant for the 
period. 

2 . Cross-Sections in Layers (2 , 4) , (4 , 6) , (6 , 8) and (8,10) 

The following definitions apply for the four layers identified 
in the present section heading and are further identified in Fig. 9. 

All of the heat transfers shown in these layers are assumed to be of 
radiative character only. The daily-averaged radiative heating (cooling) 
rate in layer (2,4) is given by 

BAL 24 = 024 - F24 (5-9) 
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Figure 10(a). 125W Longitudinal cross-section, tropical section. Refer to Fig. 9 for key. Values 

computed from data for 16 January 1974. 
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10(b)'. 12 5W Longitudinal cross-section, higher latitude section. Refer to Fig. 9 for key 

computed from data for 16 January 1974. 
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11(b). 170W longitudinal cross-section, higher latitude section. Refer to Fig. 9 for key 

computed from data for 16 January 1974. 
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Figure 12(a). 145E Longitudinal cross-section, tropical section. Refer to Fig. 9 for key. Values 

computed from data for 16 January 1974. 
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12(b). 145E Longitudinal cross-section, higher latitude section. Refer to Fig. 9 for key. 

computed from data for 16 January 1974. 
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Figure 13(a). 35W Longitudinal cross-section, tropical section. Refer to Fig. 9 for key. Values 

computed from data for 16 January 1974. 
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13(b) 0 35W Longitudinal cross-section, higher latitude section. Refer to Fig. 9 for key 

computed from data for 16 January 1974. 
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Figure 14(a). 125W Longitudinal cross-section, tropical section. Refer to Fig. 9 for key. Values 

computed from data for 16 July 1974. 
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14(b). 125W Longitudinal cross-section, higher latitude section. Refer to Fig. 9 for key 

computed from data for 16 July 1974. 
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15(b). 170W longitudinal cross-section, higher latitude section. Refer to Fig. 9 for key 

computed from data for 16 July 1974. 
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Figure 16(a). 145E Longitudinal cross-section, tropical section. Refer to Fig. 9 for key. Values 

computed from data for 16 July 1974. 
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16(b). 145E Longitudinal cross-section, higher latitude section. Refer to Fig. 9 for key 

computed from data for 16 July 1974. 
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Figure 17(a). 35W Longitudinal cross-section, tropical section. Refer to Fig. 9 for key. Values 

computed from data for 16 July 1974. 
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17(b). 35W Longitudinal cross-section, higher latitude section. Refer to Fig. 9 for kev 

computed from data for 16 July 1974. 



where 



Q 24 = daily-averaged solar absorption in layer (2,4) 

and is defined relative to A24(t) by a cosine transformation similar 
to Eq. ( 5—3 ) , and 

F24 = terrestrial cooling rate in (2,4) . 

Similarly the daily-averaged radiative heating (cooling) rates in 
layers (4,6), (6,8) and (810) are given by 



BAL46 = Q46 - 


F46 


(5-10) 


BAL68 = Q68 - 


F68 


(5-11) 


BAL810= Q810- 


F810 


(5-12) 



respectively . 

3 . Cross-Section at Air-Sea Interface (k=10) 

The radiative balance at the earth* s surface (BALB) is as 
defined in the following equation: 

•k 

BALB = QA3G - F . (5-13) 

QABG is the 24-hour average solar insolation absorbed by the surface 
as follows: 

QABG = QABG (t) (Cos z/Cos z) . (5-14) 

D. MERIDIONAL CROSS-SECTIONS OF THE VERTICAL RADIATION BALANCE 
Figs. 10, 11,... 16, 17 as previously explained represent the 
single time step of heating computations for each of the four meridians 
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used in this study. The eight figures have been divided into (a) tropi- 
cal results and (b) mid-to-high-latitude results for the mid-winter and 
mid-summer cases respectively. While the results depicted in these 
cross-sections are exhibited as representing daily-averaged values, they 
are actually based upon radiation-computations at the specific map times 
of OOOOGMT and 1200GMT on 16 January and 16 July 1974. Therefore, for 
these results to be meaningful as a stepwise part of the heat package 
subroutine of FNWC, the solar radiative absorption and reflectance terms 
would have to be recoverable as a function of GMT, i.e., 

F ( 2 , t) = QAVE* (Cos z/Cos z) (5-15) 

REF ( t) = QREF* (Cos z/Cos z) (5-16) 

etc. Thus solar disposition terms may then be utilized in connection 
with the one-hour stepwise application of the thermodynamic equation of 
the set of primitive equations used in the FNWC prediction process, 
assuming the "2/3-CL" parameterization of Eq. (2-2) . 

E. COMPARISON OF THE FOUR- LAYER WITH THE TWO-LAYER FLUX DIVERGENCE MODEL 
In previous studies (A, B, C, D) which considered only the two-layer 
flux-divergence model [layers (2,6) and (6,10)] it was assumed that the 
radiative-cooling rate in layer (2,6), BAL26, was uniformly distributed 
over layers (2,4) and (4,6). A similar remark is applicable to BAL610 
relative to the sublayers (6,8) and (8,10). In an effort to determine 
the vertical resolution of the atmospheric flux-divergences separately, 
the four values of BAL24 , BAL46 , BAL68 and BAL810 have been averaged 
over all gridpoints for the four computational dates. 
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The computed values of BAL24 and BAL68 were then compared with one- 



half of BAL26 and BAL610 respectively to determine the relative per- 
centage of heating (cooling) in these layers, as depicted in Table X 
for the mid-seasonal winter and summer dates. It is noted that the 
values of BAL46 and BAL810 are identical in magnitude but of opposite 
algebraic sign to BAL24 and BAL68 respectively. 

The results of these four-layer computations (Table X) indicates 
that in layer (2,4) there is approximately thirty percent less cooling 
than would have been deduced in the simpler two-layer model. This 
implies less radiationally induced instability in the upper troposphere 
than that originally estimated. In the layer (4,6) there is compensat- 
ing increased radiative de-stabilization. This radiative difference 
in the layer (4,6) appears to be primarily a result of the placement of 
clouds in both models. For instance, the upper clouds have tops at 
k=4, and IR net flux acts as the chief source of increased de-stabili- 
zation in layer (4,6) (and of decreased de-stabilization in layer (2,4)). 
Likewise the lower clouds have tops at k=8; so, it follows that IR net 
flux acts as the chief source of increased de-stabilization in layer 
(8,10), and finally, of reduced de-stabilization in layer (6,8). 

The tuning of solar reflectances has not introduced a significant 
radiative influence on the results of the de-stabilization differences 
just outlined. Rather it was the placement of clouds relative to 
layer boundaries that affected primarily the IR net-flux calculations. 
The results summarized in Table X are then to be considered valid if 
the cloud-parameterization model approaches climatological reality. 
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WINTER RESULTS SUMMER RESULTS 

2-Layer 4-Layer 2-Layer 4-Laye: 

Layers model model model model 
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vi. the zonal distribution of radiational 

BALANCE TERMS OF THE OCEAN-ATMOSPHERE SYSTEM 



A. GENERAL INTRODUCTION: ZONAL CROSS-SECTIONS 

The zonally distributed cross-sections of the radiation contributions 
over the ocean-troposphere system are presented in Figs. 18, 19, 20 and 
21 for winter, spring, summer and autumn seasons, respectively. The 
cross-sections, show the results after averaging over the four meridians 
on the mid-seasonal dates considered in this study. The results are dis- 
played in the format of Fig. 9. 

In obtaining zonally distributed seasonal means of radiative heating 
rates, denoted Q (<f> ) , at each five-degree multiple of latitude <j) in the 
range 20S,... 65N, values of each radiative parameter listed in Fig. 9 
were arithmetically averaged over the four meridians. At <f> = 65N there 
was only one contribution to the zonal average, while in the Southern Hemi- 
sphere latitudes 20S, 15S and 10S, only two values (on X = 125W and X = 

35W) of each parameter contributed to the means. At 5S, there were three 
meridional sets of radiative parameters entering into Q(<f>). Otherwise, 
there were four seasonal values of each radiative parameter entering into 
the arithmetic mean Q (40 at each latitude (J> of the cross-sections, Figs. 

18 , • . . , 21 . 

Consequently, near the northern and southern boundaries of Figs. 18, 

. ..,21, the listed seasonal values may not be as representative as in 
mid-latitudes. Nevertheless, the general equatorial- to-polar trend in 
the radiative-change terms appears reliable in both hemispheres. 
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18(a). Zonally-averaged radiational cross-section for tropical latitudes. Refer to Fig. 9 
. All values listed are daily averages in ly min - -*- and are computed from data for 16 January 
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Figure 18(b). Zonal ly-averaged radiational cross-section for higher latitudes. Refer to Fig. 9 for 
key. All values listed are daily averages in ly min -1 and are from data for 16 January 1974. 
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re 19(a) „ Zonally-averaged radiational cross-section for tropical latitudes. Refer to Fig. 9 
key. All values listed are daily averages in ly min"-*- and are computed from data for 16 April 
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20(a). Zonally-averaged radiational cross-section for tropical latitudes* Refer to Fig. 

. All values listed are daily averages in ly min“l and are computed from data for 16 July 
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Figure 20(b). Zonally-averaged radiational cross-section for higher latitudes. Refer to Fig. 9 for 
key. All values listed are daily averages in ly min -1 and are from data for 16 July 1974. 
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21(a). Zonally-averaged radiational cross-section for tropical latitudes. Refer to Fig. 9 
• All values listed are daily averages in ly min -1 and are computed from data for 16 October 
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Figure 21(b). Zonally-averaged radiational cross-section for higher latitudes. Refer to Fig. 9 for 
key. All values listed are daily averages in ly min -1 and are from data for 16 October 1973. 



A zonal cross-section depicting the mean annual radiative distribu- 
tion has been constructed, Fig. 22, by arithmetic averaging over the 
four mid-seasonal sets of results (Figs 0 18,.. .,21). 

It has also been convenient in earlier sections (cf. Tables III, IV, 
VIII) to compute "weighted averages" with respect to latitudes in the 
Northern Hemisphere. Thus a definition of a weighted-average parameter 
which takes into account the number of observations k. available at each 

l 

latitude has been given (Meyers, 1975) as 



14 4 




Here Q.. is the Q-value on meridian j at latitude 6. and k. = 1,...,4 is 
x ji ii 

the number of meridional observations available for the arithmetic 
average Q((Jk). Note that i = 1,...,14 corresponds to the 14 latitudes, 
(j)^ = 0, . . . , 65N . 

B. ANNUAL RADIATIVE BALANCE FOR THE EARTH-TROPOSPHERE SYSTEM 

For purposes of this summary, the term BALT (at k=2) in Fig. 9 has 
been redesignated for simplicity as R^_ (cf> ) , and annual values of R^ have 
been plotted in Fig. 23 as a function of latitude. R represents the 
net radiative- transfer rate across the tropopause, that is at the top of 
the troposphere-ocean system considered here. Similarly, the net flux 
at the surface previously denoted BALB in the key, Fig. 9, is 
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ire 22(a)* Zonal-annual radiational cross-section for tropical latitudes. Refer to Fig. 9 
key. All values listed are daily averages in ly min~^ and are computed by arithmetic averaging 
■ the four midseasonal results (Figs. 18 ,...,21). 
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e 22(b)* Zonal-annual radiational cross-section for higher latitudes* Refer to Fig* 9 for 
All values listed are daily averages in ly min - ^ and are computed by arithmetic averaging 
the four midseasonal results (Figs. 18, ***,21). 



redesignated here following Malkus (1962) as R(<f>) and its annual distri- 
bution has been graphed against latitude in Fig, 23. Finally, the rela- 
tionship between and R in any column is given by 

R = R + R (6-2) 

t a 

where R^ is the overall net cooling rate or flux-divergence of the tropo- 
spheric column, e.g.. 



R = BAL24 + BAL46 + BAL68 + BAL810 (6-3) 

a 

In Fig. 23, values of R ($) have also been plotted against latitude; how- 

a 

ever, a simpler method of obtaining R from (6-2) , namely, 

a 

R = R - R (6-4) 

Ci Lr 

has been employed. The resulting zonal annual distributions of each of 

the three parameters, R , R and R , are shown as functions of latitude 

t a 

in Fig. 23, where they are superimposed against similar functions drawn 
from Sellers (1965) « 

The comparison which follows focuses on the annual distribution in 
the Northern Hemisphere only. It should also be noted that Sellers 1 
radiative parameters were taken from climatology at the top of the mean 
zonal atmosphere , whereas this radiative model gives corresponding results 
for level k = 2, based upon tropospheric soundings over an oceanic surface 
only. 

Throughout the latitude range (0-65N) / the radiative net flux R at 
the ocean surface is substantially greater than the climatological amount 
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Figure 23. Radiative net fluxes at the tropopause (R ) at the ocean 

surface (R) ; and net-flux divergence for the tropospheric column (R_) 

from the zonal annual results (Fig. 22(a,b)). Solid lines indicate' 

values of R. R and R after Sellers (1965). 
t a 
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of Sellers. This is attributable to the reduced model-albedo values of 

Sec. IV. F. This in turn allowed more insolational-absorption in the 

ocean thereby increasing R over all latitudes as compared to Sellers 1 

values (Fig. 23). The term R-R (which equals -R ) may be regarded as 

t a 

an atmospheric "greenhouse effect" in contributing to warming of the sur- 
face both for the model and for Sellers 1 climatology. For each latitude. 
Fig. 23 indicates that the "model greenhouse" is in close agreement with 
that of Sellers. 

The Redistributions of both Sellers and of the present model are in 
good agreement across the entire range of Fig. 23, both representing cool- 
ing rates in the troposphere for all <J>* A minimum value of R^ at latitude 
20N appears on both curves and is presumably due to a local maximum of 
IR net-flux divergence F^* - (resulting from the minimum opaque 

cloud cover ) in subtropical latitudes. From SON to 65N, an increase in 
the model values of R ^ appears in Fig. 23 and seems to indicate a reduced 
IR net-flux loss associated with high CL-values (Table III) at these 
northerly latitudes. Sellers* results indicate the opposite trend, i.e., 

increasingly negative R -values from 50-60N, but his climatology may not 

a 

include the kind of detail that would indicate reasonable oceanic cloudi- 
ness over these latitudes. 

The radiative model gives R^-values which are generally larger, both 
over tropical latitudes as well as over the northern latitudes, than the 
R^-values of Sellers. As in the R((f)) comparison this is a result primar- 
ily of the reduced model-albedo. In the mid-latitude zone 25-40N of 
Fig. 23, close agreement occurs between the two R^-curves presumably due 
to such similar climatological effects as for example, the mean polar 
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front. The positive- to-negative R crossover-point occurs very near to 
40N for both Redistributions of Fig. 23. The model weighted-average is 
R^_ (cj)) = 0.048 ly min ^ if averaging is considered only to 65N. If nega- 
tive R annual values exist in latitudes (p 65N, an ocean- troposphere 
radiative balance should very nearly be established over the Northern 
Hemisphere. However data for the model were not analyzed over ice-covered 
regions poleward of 65N, where presumably higher surface albedos and 
cooler drier soundings would enable R to assume larger negative values. 

The radiative-model values of R ((})) , R (4> ) and R (c|) ) are in agreement 

t a 

with present albedo and IR net-flux observations from recent satellite 
climatology. It should be noted that Sellers 1 radiative results are 
based on climatology of London (1957) , and include also data from the 
land areas of the Northern Hemisphere; on the other hand the radiative 
model values have been based only upon soundings over the ocean. 

C. CROSS-SEASONAL EFFECTS IN THE NORTHERN HEMISPHERE 
1 . On the Net Flux Across the Tropopause, R ^ 

The seasonal distribution of the tropospheric net flux is pre- 
sented as a function of latitude in Table XI (a) where the listed R^-values 
correspond to the tropopause, level k = 2. One may discern two geographic 
regimes of annual waves in Table XI (a); a polar zone (25-65N) and a 
tropical-subtropical zone (0-25N) . 

The cross-seasonal variation of R^ in the northern regime has 
the following property: R^-values are a minimum in winter increasing to 
a maximum in summer followed by a decrease in autumn. In the tropical 
zone (0-10N) , R^-values have positive peaks in spring and autumn. The 
summer minimum of R in 0-10N is presumably attributable to increased 
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cloudiness associated with the ITCZ in this latitude zone. The subtropi- 
cal region 15-25N does not have a minimum in midsummer, and may be simply 
classified as intermediate in behavior between that exhibited by R in 
the polar and the tropical regimes. These geographic cross-sectional 
effects (Table XI (a) ) tend to confirm the general validity of the radia- 
tive model. 

2. On the Sea-surface Model Balance, R 

The seasonal distribution of the net radiative flux R at the 
earth's surface is listed as a function of latitude in Table XI (b) . Again 
a polar zone (25-65N) and a tropical-subtropical zone (0-25N) are defined 
for discussing cross-seasonal effects on R. 

The model results for R are analogous to those just specified for 
R , namely: 

(i) in the polar latitudes there is a single sine wave with 
annual periodicity and maximum in midsummer; 

(ii) in tropical latitudes 0-10N, the double peaked character 
with maximum in spring and autumn appears; 

(iii) there exists an intermediate subtropical zone of R-values 
which no longer exhibits a minimum in summer. 



D. TOP OF THE ATMOSPHERE COMPARISON OF MODEL WITH SATELLITE DATA 
1 . Net Radiative Transfer Rate at Top of the Atmosphere 

The three parameters that describe the radiative flux at level 
k = 0 are associated in the following way: 



\ ~ F T 



(6-5) 



Q 



N' 



the daily averaged incoming solar flux. 



can be related to QAVE and 
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QREF (Fig. 9) by the relationship 



Q n = Q^E (1-ALB) = QAVE - QREF 



( 6 - 6 ) 



where QAVE and OREF have already been evaluated according to Eqs . (5-3) 

and (5-6) respectively and the factor 2.00/1.92 has the effect of deriv- 
ing the mean daily solar insolation at level k = 0 before ozone-oxygen 
absorption. F , the IR net flux to space, was previously discussed in 
Sec. III.D. [Eq. (3-11)], and R^, the difference between Q^ and F^, 
represents the daily-averaged radiative net flux at the top of the 
atmosphere. In order to more easily distinguish between model-values 
of R and the same parameter from satellite climatology of Raschke et al. 
(1973), R^ is further identified as R^MOD or R^RAS, respectively. QAVE 
is assumed identical both for the model and for the treatment of the 
Raschke data in this study. 

2. Zonally-Averaged Computations of R ^ 

With the use of Eqs. (6-6) and (6-5), values of Q^, F^ and R^ 
for each of the four mid-seasonal dates and by five-degree latitude 
intervals have been computed both for the model-parameters and the satel- 
lite-observed (RAS) parameters. The results of the computations both 
seasonally and for the annual mean case are presented in Table XII. 

A comparison of R^MOD with R^RAS from the zonal-annual results 
at the bottom of Table XII shows that the crossover from positive-to- 
negative values of R^ occurs near 40N for both model- and Raschke-results . 

A statistical accuracy check of R^MOD as against R^RAS data is 
conducted using the standard deviations of the difference parameter 
(R^MOD - R^RAS) over all Northern Hemisphere latitudes. A summary of 
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the mean and standard deviations of this parameter over the latitude 
range of these five cross-sections (Table XII) follows: 





16 Jan. 


16 Apr. 


16 Jul. 


16 Oct. 


Annual 


R^MOD-R^RAS 


.004 


.024 


.019 


.014 


.015 ly min 


Std. Dev. 


.028 


.037 


.053 


.039 


.039 ly min 



The zonally-averaged mean (R^MOD-R^RAS) is positive for each 
season as well as for the annual distribution (whose mean is 0.015 ly 
min . This small positive difference was anticipated from considera- 
tion of Sec. IV. F. (Table IX) where the middle- to-high latitude albedos 
were tuned slightly too small. 

The overall standard deviation in (R^MOD-R^RAS) of 0.039 ly min ^ 
is of the same order of magnitude as that of Q^MOD-Q^RAS . This compara- 
tively small difference in the case of model versus satellite climatology 
was a considerable improvement over (A,B,C,D) where untuned albedos were 
used. 



E. STRATOSPHERIC MODEL RADIATIVE BALANCE 

It is clear that the flux-convergence in the stratosphere is given 
by 



where R is given by Eq. (6-5). Note that R is identical to BALT of 
N t 

Fig. 9 which reduces to 

R = OAVE - QREF - F 2 * . (6-7) 
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Forming the difference, R - R . leads to 

N t 

R - R = .04 CAVE - (F - F *) . (6-8) 

N t ~ T 2 

The stratospheric flux-convergence, R -R^, is displayed in zonal- 
annual format in Table XIII. Based on the Northern Hemisphere weighted 
results (bottom line) , it appears that a positive net-flux convergence 
exists between 0-65N. However the trend in R^-R^ at higher latitudes 
is towards increasingly negative radiative values in the stratosphere 
poleward of 65N such that an annual mean radiative balance may be in- 
ferred. For instance, assuming that in the latitude range 65-90N, an 
average zonal-annual value of R^-R^ = -0.0502 ly min ^ exists at 77. 5N, 
then, the cosine-weighted mean is -0.0126 ly min ^ in the region pole- 
ward of 65N. Thus the summed Northern Hemisphere weighted mean would be 
zero and a radiative balance would exist in the stratosphere. 

F. ZONAL- ANNUAL NORTHERN HEMISPHERE HEAT BUDGET 
1 . Tropospheric Radiation Budget 

The radiative heating rate of a tropospheric column may be ex- 
pressed as a function of R (<$>) by the right side of Eq. (6-9) 

a 

Q + S = R + (E + H r ) (6-9) 

va a a i 

is the storage heating rate of the column and Q is the required 
flux-convergence of sensible and latent heat energies compatible with 
the heat balance at latitude cf>. R^ is the tropospheric radiative net 
cooling rate and is depicted in Fig. 23. (E + Hp) includes the latent 
and sensible heat parameters respectively, but these parameters are 
not part of this study and will be considered to be contained in Q 
of Eq. (6-9) for simplicity. 



Ill 



Lat. 


R 

N 


R t 


\ " 1 


2 OS 


.1179 


.1030 


.0149 


15 


.1260 


.1085 


.0175 


10 


.1293 


.1074 


.0219 


5 


.1456 


.1248 


.0208 


0 


.1443 


.1215 


.0228 


5 


.1337 


.1097 


.0240 


10 


.1230 


.0969 


.0261 


15 


.1084 


.0826 


.0258 


20 


.0907 


.0644 


.0236 


25 


.0603 


.0367 


.0236 


30 


.0441 


.0293 


.0148 


35 


.0221 


.0113 


.0108 


40 


.0081 


.0052 


o 0029 


45 


-.0137 


-.0107 


-.0030 


50 


-.0273 


-.0190 


-.0083 


55 


-.0426 


-.0168 


-.0258 


60 


-.0566 


-.0290 


-.0256 


65N 


-.0563 


-.0106 


-.0457 


Wt. Avg. 


.0602 


.0476 


.0126 



TABLE XIII. Computed radiative net fluxes at the top of the 
atmosphere (R^) at the tropopause (R ) ; and the 
stratospheric absorption (R -R^) • All parameters 
are listed in the zonal-annual format in ly min" 1 . 
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The Northern Hemisphere cosine-weighted turns out to be 

R = - 0.1551 ly min ^ 
a 

and if it is temporarily attributable to mean storage-cooling S of the 

a 

troposphere alone, the mean storage rate corresponds to a temperature- 
change rate given by 

6t 1440 <V 

<6t } = 4.1 °C/ day (6-10) 

mb 

with AP^ = 800 nib in the troposphere. The resultant cooling rate (6-10) 

over the tropospheric depth, 800 mb, with zero lateral flux divergence, 

2 

is approximately -1.14°C per day averaged over the mean cm tropospheric 
column. 

Thus the general circulation of the atmosphere would have to 

function to bring about the atmospheric heat balance. This process 

would require the proper flux-convergence Q to offset the annual radia- 

2 

tive loss of the average cm column of the troposphere. 

2 . Zonal-annual Heat Budget of the Ocean 

An equation similar to (6-9) may be written for the zonal-annual 
heat budget for the ocean as a function of R(cj)) . 

Q + S = R - (E + H r ) (6-11) 

vo o r 

Here, S is the storage rate of the ocean water mass and Q is the 

O vO 

required mean oceanic heat flux divergence for a balance at latitude <J). 

The primary ocean mass heating parameter is the net radiative heating 

flux at the ocean surface, R (Fig. 23), and (E + Hp) is for convenience 

considered part of O 

*vo 
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The Northern Hemisphere annual weighted mean value of R is 
R = 0.2027 ly min 1 

which corresponds (Table XI) to a mean radiative heating rate in the 

water-mass column. It should be noted that this mean heating rate is 

noticeably larger than R and of opposite sign to the corresponding 

a 

tropospheric radiative cooling effect (R = -0.1551 ly min 0 How- 

a 

ever, if the three-dimensional oceanic flux divergence of both 
sensible and latent heat is considered, the oceanic heat balance should 
result. 



114 



VII. CONCLUSIONS 



The major change in this radiative model from the model evolving 
from (A,B,C,D) was the introduction of tuned cloud reflectances. This 
innovation brought the resultant model-albedos into agreement with 
satellite climatology for comparable data dates. 

Gridpoint comparisons at the top of the atmosphere of model net flux 
and of the net incoming model insolation were made with similar 
parameters from Raschkes 1 satellite climatology (1973). Closer agree- 
ment was obtained between Q MOD and Q RAS than in any of the preceding 

N N 

studies (A,B,C,D) . Similarly, close agreement in a least-squares 
sense was achieved here between F MOD and F^RAS primarily because of 
an improved formulation for the F - computation than was used in 
earlier studies. Finally the ocean-atmosphere system net flux 

\ “ 0 N - F t 

has only a small least-square error when model and satellite results 
are compared with properly stratified gridpoint and time cross-sections. 

In this study, mid-latitude cloud reflectances were also tuned to 
obtain closer agreement between the global albedo and corresponding 
values from satellite climatology. The result was to systematically in- 
crease relative to Raschke's values. This result was undoubtedly 
due to a mismatch between the model clouds and those involved in the 
climatological average. The conclusion reached is that if tuning of 
cloud reflectances of this nature is to be attempted in mid-latitudes 
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using FNWC soundings, more accurate cloud-parameterization formulas 
should be established against observed cloud-cover amounts for identical 
experiment-periods similar to GATE . If then further cloud-reflectance 
tuning is necessary, it could proceed from a more certain knowledge of 
the "ground truth" provided by these experiments. 
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APPENDIX A. COMPUTER PROGRAM 



//FORT .SYS IN 
COMMON 



502 

500 



C 

r 

c 

c 

c 



CD * 

6T( 11 ) ,DA74( 11 , 5) ,CL (2 ) 

CC MMO N/ ABS / FA2 » AL PH AG * C Z , 5 ECZ , F S , FA D J 
COMMON/ 7 k I G/ SI, C7,SD,CD, AL AT 
COMMON/ VP/ Q( 11 ) , QS (11 ) 

COMMON/ ARM/ RM 
COMMON/ R ASK C/RAS VAL 

CO MMO N/ CASK/ AL BRAS (93) ,F2RAS( 53 ) , ISO 
COMMON/ FACT/ K,IS£A 
FQRMAT(4X,2(F8.4,2X> ) 

FORMAT (2X,F6.1,lX,F7.2,lX,5(F9.5,iX),<>X,2I2) 
DO 63 ISEA=1,4 

■Q.l) DEC=-21. 07917/57. 29578 
' “ DEC = 8.43333/ 57.29578 

DE C = 21.5000/57.29578 

DE C = -8.225/57.29576 

RM = . 9837 1 5**2 
RM = 1 • 30 33 3**2 
R M = 1 . C 1 644**2 
IF (ISEA. EG. 4 > RM = . 997165**2 



I F ( I SE A . 
IF ( IS EA. 
IF (ISErt. 
IF ( ISEA. 
IF (IS EA . 
IF ( ISEA. 
I F ( IScA 



2Q .2 ) 
Eg. 3 ) 
EQ.4) 
Ey.l ) 
6Q.2 5 
EQ.3 ) 



THE FOLLOWING READ : 
CLIMATOLOGY DATA C p 
ROUTINE REFT. 



IT AT EM EN T 
RASCHKE 



= 



APPLIES TO THE ... 
ET AL. AS USED IN 



TELL I 7 F 
7F!c S U B- 



1 



1C 



CC 9 I $0 = 1 » 92 

RE AD ( 5, 502 1 ALBR A S ( I SO ) , F2R AS ( I SO ) 

CONTINUE 
DC 66 150=1,93 
CO 10 1=1,11 

FEAD(17 ,500, END=99 ) ( DATA ( I , K ) , K=l, 2), 0(1), QS (I) 

(CA7A(I ,KK) ,KK=3, 5) ,LAB1,LAB2 

CONTINUE 

IF (ISEA.EQ.3.AND.ISC.EQ.25) LABI =24 
CALL ANGLE (CZ, SrC Z , I 30, LASl , DEC > 

1F( ISEA.EQ.l .AND. ALAT.LE.25. ) R AS VAL 
IF (ISEA.EQ.l . AND. A L AT “ 

. H ' ’< l ) . rA L A j | 

. Ar» D . <^ LA i i 
AND. ALA I 
tQ .3 • AND • ALa 1 
EQ.4 . AND. ALAT , 

EQ.4. AND, ALAT, 



, E Q . 2 
. EQ.2 
> c Q * 



LE, 
G i c 
LE, 
GT 



i 'j 

L 

. GT 



25. ) 
25 . ) 
25. ) 
,25.) 
25. ) 
25. ) 



I F (ISEA 
IF (ISEA 
IF ( ISEA 
I F ( I S t A . 

IF (ISEA. 

IF (ISEA. 

C -■ LL BBB 

ALPFAG= ( ,54*( ,7-CZ) )+.06 
I F ( ALPHAG.LT..06 ) ALC>K4G=.06 
FS=.651*CZ*1.92/RM 
F ADJ= (2.00*CZ )/RH 
F A 2= ( .349*1. 92*CZ) /RM 
CO 20 1=1,4 
CL(1) =0.0 
CL ( 2 ) =0 .0 

IF ( I. EQ.2. OR. I .EQ.4) CL(1)=1.000 
IF(I.GE.3 ) C L ( 2 ) = 1 . 00 0 



R AS VAL 



R A 3 V A L 
R AS VAL 
R AS VAL 
RAS VAL 
R AS VAL 
RASVAL 



, H 
,8 
i 3 5 
,65 
, 6 
, 7 
, 4 
,6 



ir(i.GE.3) uin = i.wu 
WT=(1.0~D4TA(6,3 ) )*(1. 
I F ( I . EQ . 2 ) W 7 = 0474(0,3 

7 r / T r r \ i ,.i T .. r , a t a / o *5 



WT=DATA(2,i 



( 1 .O-DAT A ( 2, 2 ) ) 

- ")*(1.Q-DATA(2,3) ) 



/ ^ l L • \J L- H I Hid, J I I 

>*(1 .G-DATA (6, 3) ) 

( ^ n / v A / a **> \ 



20 

66 

68 

95 



•D 
p WT 



T A ( 6 , 3 ) 

) 



IF (I. EQ.2) 

IF (I. EQ.4) WT = DAT A ( 2 , 3 ) 

CALL A6S0RB( CL ( 1 ) ,CL(2 ) 

CONTINUE 
CONTINUE 
CONTINUE 
STOP 
END 

SUBROUTINE A3S0RB ( Cl , C2 , WT ) 

CCMtMCN/ABS/ FA2, AL PHAG, C Z , S ECZ , FS , FAD J 

C C- M MON / R A S KC / R A S V A L 

IF(C1 .GT. .5.0R.C2.GT..5) GO 7C 20 

AL PHAR= ( ,085 + ( . 2 5 074* ALOG 10 ( SEC Z ) ) ) 

IF (ALPHAR.GT. 1.000) ALPHAR= 1.000 
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A = ( ( 1 • - A L PHAR ) * ( 1 .-ALPHA G) ) / d.- < ALPHAR* ALPFAG ) ) 
SOL=-A*F S 
GC TO 21 

2C R1=.54*RASVAL 
R2=.66*RASVAL 
IP(C1.LT..5) R 1 = 0 • 0 
I F ( C2 . L7. • 5) R2-G.G 
A = ( 1.-R1 )v(l .-R2 )*< l.-ALPHAG) 

C=1.0--( (Rl-"‘ c 2) + (R2*ALPHAG)+(Ri*ALPHAG) ) 

C=2.0*Rl*R2*ALPhAG 

SCL=<A*FS ) / ( D+C ) 

21 !F(Cl.LT..l.AND,C2.LT,.l) CALL AEOO< WT, SCL ) 

I F < Cl.GT. . 9.ANQ.C2. LT. . 1 ) CALL A61Q(WT,SCL) 

IF (Cl.LT.. 1.AM0.C2.GT..9) CALL A BO 1 ( W7 , SCI. ) 
JF(C1.GT..9.AND.C2.GT..9J CALL AB11(WT,SCL) 
RETURN 
END 

SUBROUTINE TD< ANS, IND,N1,N2 ,SEC, C) 

CCNMON BT( 11 ) ,DATA( 11, 5),CL(2i 
U = DATA(M,4)- DAT A ( N2 » 4 ) 

A =.203 

ANS=C*( (U*SEC )**A ) 

I F ( IND.GT.10. ) AN S = 1 . O-ANS 

RETURN 

Ei\C 

5 L GROUT INF BBB 

COMMON BT( 11 ) , DATA! 11, 5 ) ,CL ( 2 ) 

£ 1 GMA = 1 .170403 
E= 1 • E-7 
DC 10 1=1,11 
A = C AT A ( I ,2 ) + 273.16 
07 ( I )=( A**4) *SIGMA*E 
10 CONTINUE 

C AT A ( 2 ? 4 ) = ( DAT A. ( 1 , 4 ) +OATA ( 3 ,4 ) ) /2 .00 
D A T A ( 2 , 5 ) = ( DAT A ( 1 , 5 ) +DATA ( 3 , 5 ) J/2.CC 
RETURN 
END 

SUBROUTINE ANGLE (CZ,SECZ,IS,L2,D) 

CCNMGN/TRI 0/ ST, CT , SO, CD, AL AT 

R=57. 29578 

h=25.Q/R 

IF(IS.LE.25) H=55.0/R 
I F ( IS .GE.26.AND. IS.LE. 50) H=1C.0/R 
A=973.752 
S= 2.0 

I F( IS. GE. 26. AND. IS.LE. 50) S=I. 0 
AK=L2+1 

IFdS.GE.68) AK=63-L2 
e={S’M(32.-AK>**2) ) 

ST = ( A-6 ) / ( A+B ) 

CT= (1.- (ST**2> )**.5 
SC=SIN(C ) 

C C =COS ( D ) 

CF=COS ( H ) 

CZ=(S7*SD)+(CT*CD*CH) 

SECZ=1.00/CZ 
TANT= ST/CT 

ALAT=( ATAN(TANT) )*57. 29578 

RETURN 

ENC 

SUBROUTINE ABOO(WT,SOL) 

CGMMON/A6S/FA2, ALPHA G,CZ, SEC Z, FS, FAD J 
CALL TD(U24,0,9, 7,SECZ,.271 ) 

CALL TD(U26,0,9,5,SECZ,.2.71) 

CALL TC(U23,0,9,3,$5CZ,.27i) 

CALL TC(U21G,0,9,1,SECZ,.271) 

CALL -TC(U10,12,9,1,SECZ,.271) 

A24=FA2*U24 
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A46=FA2*< U26-U24 ) 

A68=FA24(U28-U26 ) 

A610=FA2*(U210~U28> 

T R ANA = F A 2*U10 
A I 1 0 = T R ANA*< 1.- ALP HAG ) 

R EFA=FA2-A24~A46-A6S-A81C~AI 10 

CALL REFT ( 1, A24, 446, A68, A81C, A I 1C, TRANA , REFA ,S0L , WT ) 

RETURN 

END 

SUBROUTINE A 810 ( WT , SOL ) 
CCMM3N/ABS/FA2,ALPHAG,CZ,SECZ, F$ , FAD J 
CCMMCN/ RASKC/R ASVAL 
R 1 = • 464 R AS VAL 
A 1=. 20 

CALL TC(TD24,12,9,7,SECZ,.271) 

CALL TD(TDo8,12, 5, 2,1 ,6667, .271 ) 

CALL TO (TD61 0,1 2, 5,1, 1.66667, .271) 

CALL TD (T 08 10,12 ,3,1, 1.6667, . 271 ) 

F4D=FA2*TD24 
F40=R 1- F4D 
F2L = F4li*TD24 
A24=FA2-F40+F4U-F2U 
A46=F4D*A1 
F6C=F4D*< 1.-R1-A1 ) 

F8C=F6C*T068 

F10D=F6D*TD610 

F10U=Fi CD4ALPHAG 

F6U=F1CU*TDS1G 

F6U=F10U*T0610 

F6CD=FfcU*Rl 

FfiCC=F6C0*TD68 

F 10DC= F6DQ4TD610 

Ad£ = F6D-FSD+P8’J~F6U+F6DD-F80D 

A £10= F8C-F10D+F10U-F8U+F8DS-F10D0 

AI10=(F10D4{ 1.- ALPHAS) )+F10DD 

TRANA=F 1CD+F10DD 

R6FA=FA2- A24-A46- A68-A8 10-A ! 10 

CALL REFT ( 2, A24, A46, A68 , A8 10 , A I 1C , TR ANA , RE F A , SOL , WT ) 

RET URN 

ENC 

SUBROUTINE A301(WT,S0L) 

CCMMCN/ABS/FA2, ALPHA G,CZ, SECZ, FS , FAD J 
COMMON/ R AS KC/R AS VAL 
R 2= . 504 R A S VA L 
A2=.30 

CALL TD(TD24.12,9, 7,SFCZt.271 ) 

CALL TC(TD26,12,9,5,SECZ,.271) 

CALL TC(TD28,i2,9,3,SECZ, .271) 

CALL TD( 7 0910,12 ,2,1, 1.6667, .271 ) 

CALL TC(TD68,12, 5, 3, SECZ, ,271) 

CALL Tu( T D48*i2,7, 3, SECZ » .271 ) 

FAD=FA24T024 

F6C = F A 24TD26 

FcC=FA2*7D23 

F EU = F 8C *R2 

FfcL=F8U4T068 

F4L=F8U*TD48 

F2l=F8U*TD28 

A24 = FA2— t-40+F4U— F2U 

A4fc=F4C- F6D+F6U-F4U 

A68 = F6C-F8D«-F8U-F6U 

F 9C = F 804 ( 1.-R2-A2) 

F1GD=F9C*TD910 
F10U=F1CD*ALPHAG 
F9L=F1GU*TD910 
F 9 C D= F 9 U * R 2 
F10CC=F9DD*TD910 

A£10=( F8D4A2 ) +F9D-F 10 D+F 10 U-F9U+ F 9DD-F1 0 DD 

TRANA=F 10D+F10DD 

AI10=(F10D*( l.-ALPHAG) )+F10DD 
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P E FA= FA 2- A 2 4- A 4 o- A 6 6~ A 310- A I 10 

CALL Rc FT ( 3 , A24 , A46 , 46 6 , A8 10, A I 10 , TRAN A , RE FA , SOL , WT ) 

RETURN 

END 

SUBROUTINE ABil(WT,SCL) 

C C M.MON / A 8 S / F A2 . A L P H A G , C Z , S EC Z f F S f F AD J 

CCMMCN/RaSKC/RASVAL 

R1=.46*RASVAL 

R2=.50*RASVAL 

A 1 = • 20 

A2=.30 

CALL TC(TD24,12 f 9»7,SSCZ,.271) 

CALL T C ( T 068 ,12, 5,3,1.6667, .271 ) 

CALL TC ( 7 0910 , 12 , 2 , 1 , 1 . 6667 , * 27 1 ) 

F4C = F A 2^“ D 24 
F4U=R1*F4D 
F2U=F4U*T024 
A24=FA2-F40+F4U-F2U 
A 4 6= F40*A 1 
FCC=F4D*( 1.-R1-A1 ) 

F8C=F6G*TD68 

F6L=F8C*R2 

F6l=F8U*TD6e 

F6CD=F6U*R1 

F8DD=F6DD*TD68 

A68=F6C-F8D+F8U-F6U+F6DC-F3DD 
F9C=F8D*( l.-R2-A2)+F8DD*< 1.-A2 ) 

F1GD=F9C*7D910 
F10U=F10D#ALPHAG 
F 9 1 = F 10 L * T D9 1 0 
F9DD=F9U*R2 
F 1CDC=F9DQ*7D91G 

A8 10= ( ( F 3 0 + F800 ) *A2 )+ F9D-F1CD + F 1GU-F9U + F9DD-F 10DD 
Ai 10= ( F!OQ*< l.-ALPHAG) )+FluDD 
TPANa=F10D+F10DD 
REFA=FA2-A24-A46-A68-A810-AI 10 

CALL kEFT(4,A24, A46,A68,A810, A I 10, TRAN A, RE FA , SOL , WT ) 

RETURN 

END 

Ste POUT INF REFT ( NSUB, A24 , A46 , A68 , A 8 10, A I 10 , TRAN A , 
1REFA ,SOL , w7 ) 

COMMON ET ( 11 ), DATA (11, 5) ,CL<2> 

COM MO N / A B S / FA2, ALPHAG,CZ, SECZ, FS, FAOJ 
COMMON/ TRIG/ ST , CT , SO , CD , AL AT 
CCMMCN/ ARM/ RM 

COMMON/ RA.SK/ ALBRAS ( 93 ) ,F2RAS(93) ,ISO 

CCMMC N/PASKC/ RASVAL 

COMMON/ FACT/ K,ISEA 

I F ( NSUB . GT • 1 ) GO TO 10 

C A 24=0 • 0 

CA46=0. 0 

C A63= 0 . 0 

CAE 10=0 • 0 

CAI 10=0.0 

CS 110=0.0 

C P EF = 0 • 0 

10 CA24=CA24+A24*NT 
C A 46 = C A 4 6 + A46*WT 
CA68=CA68+A6S*WT 
CA£10=CA810+A810*WT 
CAI 10=CA1 10+AI10*WT 
CSI 10=CSI 10+SDL*WT 
CR EF = CR E F+ ( ( FS-SOL ) *w T ) +R EF A*WT 

CALL IR (NSUB, A10, A8, A6, A4,A2, FID, F8, F6, F4, F2, WT, FT ) 

I F ( NSUB • L T . 4 ) RETURN 
AEG=(CAI10+CSI10J/FADJ 
T AND= SC/CD 
TANT=ST/CT 
T AN = -TAND*TANT 
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H=ARCOS (TAN) 

S l NH= S I N ( H ) 

BARCCS= ( ( ST ^-SD# H ) + (C D*CT *S INH ) )/3. 14159 
F 200= 1 . 92/ RM 
0 A VE= F200*3ARC CS 
A F 2 =- F 2 

CF 5F= (CR5F*QAVE) / ( FAD J*. 9600 ) 

B ALT = CAVE-QREF-F2 

0 24= < CA 24*QA VE ) / ( F AD J*. 9600 ) 

A F24= -( F2-F4 ) 

BAL24=C24+AF24 

Q46= ( CA 46*Q AVE ) / ( FAD J*. 9600 ) 

A F 4 6= - ( F4-F6 ) 

BAL46=G46+AF46 

06 6 = ( CA 6 8* 3A VE ) /{ FAD J*.96C0 ) 

A F6 S = - ( F6-F8 ) 

B A L 6 8= Q68 + AF68 

0810= (CA81C*QAVS >/( FAD J*. 9600) 

AF810=-(F8-F10) 

8 AL 810= 08 10 + AF810 
CA 8G= { ABG *QAVE) /.9600 
B A L 6= CABG-Flu 
AF10=-F10 
CA26=CA24+CA46 

026 = ( ,5*CA26*GAVE) /( FAD J*. 9600 ) 

A F 26=- . 5* ( F2-F6) 

GAl26=C2o+AF26 

CA61C=CA68+CAB10 

Q6 10= ( . 5*CA610*QAVE)/(FADJ*.9600 ) 

AF610=-.5*(F6-F1Q) 

6AL610=C61C+AF610 

CIF24=B^L24-8AL2o 

DI F46=BAL46-SAL26 

DI F66 = BAL68-BAL610 

CIF810=BAL310-EAL610 

AL E = CRt F/FAD J 

ALBDI F= A LB-ALB RAS ( ISO) 

F 2 C I F = FT - F 2R AS (ISO) 

RNRA$=(CAVE*1»G41667*(1.-ALBFA$ ( I SO ) ) ) -F2R AS ( IS 0 ) 

R N i v 0 C = < CAVE*1. 041667* ( 1.-AL0) )-FT 

D I FRN = P NMCD-RNRAS 

RETURN 

END 

SUBROUTINE IR(NSUB,FlG,F8,F6,F4,F2,SlG,c8,E6,£4,E2, 
1UT ,ET ) 

IF(NSUB.GT.l » GO TO 10 
OF F = 1 .0 
810=0.0 
EE=C.O 
E 6=0 .0 
E4=0.0 
1 2=0 • 0 
ST =0*0 

10 CALL FF101 Flo, OFF ) 

CALL FF 8 ( F8 ? OF F I 
CALL F F 6 ( F 6 » 0 F F ) 

CALL FF 4 { F 4 » OF F ) 

CALL FF2( F2» JFF ) 

CALL FFTOP ( FT » OF F ) 

E10=E10+WT*F10 

E8=68+WT*F9 

E6=E6+wT*F6 

E4 = E4tWT *F4 

E2=E2+WT*F2 

E7=ET+WT*FT 

1 F (NSUb.LT. 4) RETURN 
E1G=E10/1 440 . 

E6=E8/144G. 

E6=E6/ 1 440 • 

E4 = E4/ 1440 • 
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E2=E2/1440. 

F T = E 7 / 1 4 4 0 • 

RETURN 

EKC 

C SUBROUTINE FF 10 { FIO, OFF ) 

COMMON BT( 11 ) , DATA( 11,5) ,CL{ 21 
I F ( CL ( 1 ).GT..5.0R.Cl(2).GT..5) C-0 TO 10 
CALL E!nC(EW8,GFF,3,1) 

CALL El-.C( EW6,GFF, 5,1) 

CALL EWC(EW4,0FF,7,1) 

CALL EKC ( EW2 ,3FF ,9,1 ) 

CALL EwCTEWl, OFF, 10,1) 

CALL EWC(cW9,0FF,2,l> 

CALL BOUND ( WAVE , 0 FF , 1 1 , 1 ) 

0 1= ( BT ( 3 ) - 37 ( 5 ) )*(EW8+EWo) 

B 2= ( BT ( 1 ) — BT ( 3 ) ) *EW3 

Al = < ( Bit 1)-BT( 5) ) - • 5 * ( 31+92) ) 

fil = BT( 10)*'wAVF 

B 2 = ( BT ( 9 ) -B 7 ( 10 ) )7(EW2+EW1) 

E 3= ( BT ( 7 ) - BT ( 9 ) ) * ( £ W4 + 5W 2 ) 

B4=(BT(5)-3T(7) )*<EW6 +EW4) 

A 2= ( BT (5 I - ( .5* (81 + 32 + 33 + 64) ) ) 

B1=(1.-(.5*EW9) ) 

A 3 = ( 67 ( 11-61(2) ) 7*81 

1C F1Q*U l.-CL< 2) )*A1 ) + ( ( l.-CL< 2) ) * (l.-CL( 1) )**2) + 

1 (CL (2J7A3 ) 

RETURN 

EKC 

r 

SUBROUTINE FF8(F3,Q P F) 

COMMON ET<11),DATA<11,5) ,CL<2) 

I F<CL(1 ).GT..5.0R.CL (2) .GT..5) GC TO 10 
CALL EhC(EW68,0FF,5,3) 

CALL EWC<EW48,CFF,7,3) 

CALL EWC(£W28,0FF,9,3) 

CALL EWC(EW18, OFF, 10,3) 

CALL EWC ( cW810,OFF, 3, 1 ) 

CALL BCl'ND(WAV£, OFF, 11,3) 

B1=EW68*(B7( 3)-3T( 5) » 

B 2 = ( F W 6 8 + E W4 8 ) * ( 37 (5) -37(7) ) 

E3=(£W48 + EW2 8)*( BT( 7) -BTC 9) ) 

B4=(EW28+EW18)7( 3T ( 9 ) - BT (10) ) 

85= WAV E73T ( 10 ) 

A 1=BT (3 )-( .57(81 +62+33 +B4+B 5) ) 

A2=(BT( 1 ) - BT C 3 ) ) *( 1.- ( . 5*EW81C) ) 

A3=(BT( 3)-BT( 5) )*( l.-< .5*EW66) ) 

10 F8=( ( 1 . -CL ( 1 ) I7A1 )+( (l.-CL(l) )*<1.-CL(2I )7A2)+(CL(1J 
17A2)+(CL( 1)*( l.-CL (2) 5*A2) 

RETURN 

EKC 

C 

SLBRCUTINE FF6(F6,0FF) 

COMMON BT C 11 ) , DATA (11, 5) ,CL(2) 
IF(CL(1).GT..5.0R.CL(2).GT..5) GO TO 10 
CALL EwC(EW68,0FF,5»3) 

CALL EaC( EW46,0FF , 7, 5 ) 

CALL EWC(EW26,0FF»9,5) 

CALL EWC(EWlo,0FF,10,5) 

CALL EWC ( tW6 10 ,0F F , 5 ♦ 1 ) 

CALL 3CUND(wAVE,0FF, 11 , 5) 

B 1 = W A V 1 7 8 T (10) 

G2=(3T( 9 ) - B V ( 10) ) * < EW26 + EW16) 

B3=(3T( 7) -BT (9 ) )7( EW46+EW26) 

34=EW467(BT(5)-BT(7) ) 

85=EW63*(67(3)-3T( 5) ) 

A1=<6T(3)-(.5*(B1+B2+B3+B4+B5) ) ) 

A2=(BT( 1)-8T(3) ) * ( 1 . - ( . 57 ( 5 W63+5W610 ) ) ) 

A 3 = ( 1 • - ( . 57EW68) ) * ( BT ( 3 ) - 3T ( 5 ) ) 

10 FC=( (l.-CL(l) ) 7 A 1 )+( (l.-CL(i) )7(1.-CL(2) )7A2) + (CL( 1) 
17A3 ) + (CL( 1 ) 7 ( l.-CL (2) )*A2) 

RETURN 
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ENC 

C 

SUBROUTINE FF4 (F4,0FF) 

COMMON BT( 11) ,DATA( 11 , 5) ,CL (2) 

IF(CL(1 ).GT..5.0R.CL(2) .GT..5) GO TO 10 
CALL EhC(F.W4o,0FF,7,5 ) 

CALL EWC( EW48,0FF , 7, 3 ) 

CALL EwC< SW24,0FF, 9,7} 

CALL EWC(EW14,CFF, 10,7) 

CALL EWC(£W410,0FF,7. 1 ) 

CALL BGUN0(WAVE,0FF,11,7) 

WAV6=WAVE*57 ( 10) 

61*(BT( 5 )-6T ( 7 ) i *EW46 
E2=(BT < 3 ) - B T ( 5 ) ) * ( EW46 + EW48 > 
fi 3= ( BT ( 7 ) - BT ( 9 ) ) *EW24 
B4=(BT (9J-BT (10) ) *( EW14+EW24) 

A1 = BT (3 )-( ,5*( 61+B2 + 83 + B4 + WAVE ) ) 

A2 = 6T ( 7 )- ( . 5 *( 83+ B 4 + WAVE) ) 

A3=(bT( l)-BT(B) J *( 1„-U5*(EW48 + EW410) I ) 

10 F4=( ( l.-CL (1 ) )*A1 ) + {CLU)*A2 ) + ( ( l.-CL( 1 ) )*( l.-CL( 2) ) 
l*A 3 ) 

RETURN 

ENC 

C 

SUBROUTINE FF2(F2,0FF) 

COMMON BT 1 11 ) ,l)A 7A( 11,5) , CL(2 ) 
IF(CL<li.GT..5.C'R.CL(2),GT..5) GO TO 10 
CALL EWC ( EW2*,GFF, 9, 7 ) 

CALL EWC ( EW26,f)FF,9, 5 ) 

CALL EWC(cW26, OFF, 9, 3) 

CALL EWC(EWi2, OFF, 10,9) 

CALL EWC( EW210 ,OFF, 9, 1 ) 

CALL 8 C U N D ( WA V t , 0 F F , 1 1 , 9 ) 

C1 = (BT( 1 > — 3T ( 3 ) ) 

C2*(l.-( ( EW28+EW21O/2.0) ) 

A3=C1*C 2 

61 = V.AVE-BT ( 10) 

62=EW12*( B T ( 9)-B7( 101 ) 

B 3 = ( EW26+EW28 ) * ( BT ( 3) -3T ( 5) ) 

B4 = ( E W2 4+ c fa 2 6 i * ( 37 (3)-BT (7 ) ) 

B 5 = ( EW24)*(6T( 7)-BT(9) ) 

B=61+ 32+83+84 +B 5 
Al=(87<3)-(.5*6) ) 

B3--EW24*! BT ( 7)-8T< 9 ) ) 

A 2 = ( 87 (7)-(.5* 131+82+33) ) ) 

10 F2=<A3*<1.-CL<2) )*<1.-CL( I) ) )+<Al*<l.-CL (im + (A2* 
1CL( 1) ) 

RETURN 

END 

C 

SUBROUTINE FFTOP ( FT , OFF ) 

COMMON BT(il) , DAT A ( 1 1 , 5 ) ,CL<2) 

IF(CL(1).G7. .5. OR. CL(2).G7..5) GO TO 10 
CALL EWC(BW01 ,JFF, 11, 10) 

CALL EWC(Efc'C2, OFF, 11,5) 

CALL EWC(EW04, OFF, 11,7) 

CALL EWC ( E wOb , OFF , 11 , 5 ) 

CALL EWC(EW03, OFF, 11,3) 

CALL EWC(EW010, OFF, 11,1) 

B1=<EW01+EW02)*< ST (9) -81 ( 10 ) ) 

B2=( EW02+5W04)*( 3T( 7) -3T(9) ) 
B3=(EW04+EW06)*(BT(5>-87(7) ) 

B4=( EW06+SWG8)*( BT(3) -37(5) ) 

A1 = (BT (3)-(.5*(Bl+32+63+B4) )) 

B3=( l.-( .5*( EW08+EWG10 ) ) ) 

A2 = 83* ( BT { 1 1 -BT ( 3 I ) 

A2= ( BT ( 7 ) - ( • 5* ( B 1 +82 ) ) ) 

10 FT=(A1* ( l.-CL( 1 ) ) ) + ( A2*<1.-CL ( 1 ) )*( l.-CL ( 2) ) ) + ( A3* 
1CL( 1) ) 

RETURN 

END 
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SUBROUTINE BOUND ( WAVE ,CFF,N1, N2 ) 
COMMON BT( lli ,DATA(11,5) ,CL(2) 

7 = CAT A ( 10, 2) +273. 16 
U = C AT A ( N 1 » A ) - 0 A f A { N 2 , 4 ) 

IFtU.LT. .00001 ) U= . J0C005 
AL=AL0G10 (U) 

D= . 353* AL-.44 
A 1 = 6. 34 V) **Q 
D=~. 03455* AL-. 705 
A2=U**D 

A3=(6.QG/( .35 3* AL + 3.56) ) 
W7VE=A1*A2~A3 
IF(OFF.GT.10.) RETURN 
C = CATA(N1, 5)-DAT A (N2» 5) 

D = ( U+.028o) **.26 
81=.07262*(l.-(. 62556*0) ) 

82 = AL0G10 (C )+( 1 .064) 
h A V E = w A V t + 3 1 *6 2 
RETURN 
END 

SUBROUTINE E WC ( E W , OFF , M 1 , M 2 ) 
COMMON 31(11), DATA ( 1 1 , 5 ) , CL ( 2 ) 

U = DAT A ( M 1 , 4 ) - 0 AT A { M2 , 4 ) 

D = ALC'G1C(U + .Q1C) 

A 1= ( • 240*D ) + . 622 
Efc=Al 

IFC3FF.GT.10.) RETURN 
C = CA7 A( Ml , 5) -DAT A( M2, 5 ) 
81=l.-(.62556*( (U+.0286)**.26) ) 
B2=AL0G1 0(0+1. 0 64 
EU=Al+( .07262*31*62 ) 

RETURN 

ENC 



CONTROL CARD FORMAT FOR 2/3 CLOUD DATA AND "RACIATIVc 
SOUNDING" DATA ON NPS-304 

//GO . FT 17F001 DC UN IT= 3400-4, VOL = $S D =N PS204 , D I $ p = (OLD, K 
// DCB = ( DEN* 2, RECFM = FB, LRECL= 77, BLKS I ZE = 84 H, LABEL* ( I,S 
// DSN=UCGSE7 
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